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I.  Introduction 

The  development  of  the  theory  of  causal  operators  defined  on  a 
resolution  space  was  initiated  a decade  ago  in  response  to  the  failure 
of  classical  Hilbert  space  methods  to  yield  a solution  to  the  quadratic 
optimization  problems  of  mathematical  system  theory.  Although  this  theory 
has  now  achieved  a considerable  degree  of  maturity  yielding  viable 
solutions  to  problems  arising  in  network  synthesis  Cl],  feedback  system 
stability  [9],  sensitivity  theory  [1],  and  stochastic  processes  [7,8], 
the  solution  to  the  original  quadratic  optimization  problems  has  remained 
elusive.  Porter  and  DeSantis  [12,23]  have  solved  a deterministic  servo- 
mechanism problem,  Stienberger,  Silverman,  and  Schumitzky  [22],  have 
solved  a deterministic  regulator  problem  and  Saeks  [7],  has  solved  a 
stoch.istic  identification  problem.  The  original  goal  of  a general  theory 
for  quadratic  optimization,  however,  has  yet  to  be  achieved. 

The  purpose  of  the  present  paper  is  to  present  a derivation  of  the 
Wiener-Hopf  filter  using  resolution  space  techniques.  Although  still  re- 
stricted, we  believe  that  the  tools  employed  are  indicative  of  the  tech- 
niques which  will  eventually  lead  to  a general  theory  of  quadratic  opti- 
mization. Indeed,  even  for  this  restricted  filter  the  deriviation  re- 
quires several  recently  developed  results  from  the  theory  of  causal 
operators.  These  include: 

i)  the  additive  decomposition  theorem  for  Hilbert-Schmidt  oper- 
ators [10], 

ii)  the  miniphase  factorization  theorm  [7,  19], 


iii)  the  theory  of  resolution  space  valued  stochastic  processes 
[7,  8,  11], 


iv)  the  quasi-nilpotence  theorem  for  strictly  causal  operators 


[9,  10]. 


Figure  1.  Wiener-Hopf  Filter 

A basic  filtering  problem  is  illustrated  in  Figure  1,  where  the  signal 
"X”  and  noise  "n"  are  mixed  together.  The  goal  is  to  pass  the  mixed  signal 
through  a filter  T (to  be  designed)  to  get  an  output  y such  that  y is  the 
"best  copy"  of  X that  can  be  achieved.  When  X and  n are  assumed  to  be 
stationary,  zero-mean  and  independent  stochastic  processes  and  the  "best 
copy"  of  X is  defined  in  the  sense  that  the  error  e (*  X - y)  has  a mini- 
mum variance,  the  solution  is  nothing  but  the  classical  Wiener-Hopf  filter. 
[2]  In  this  paper,  we  would  like  to  formulate  this  problem  in  Hilbert 
resolution  spaces,  i.e.  we  assume  X and  n to  be  Hilbert  space  valued  random 
variables  (zero-mean  and  independent)  and  attempt  to  find  a filter  T 
(causal  operator)  such  that  the  output  y is  the  best  copy  of  X in  the 
sense  that  the  "variance  operator"  of  the  error  is  minimal. 

Before  we  continue,  let  us  explain  the  terminology  we  have  just  used. 
Hilbert  Resolution  Space 

By  a Hilbert  Resolution  Space,  we  mean  a 2-tuple,  (H,E),  where  H is 
a Hilbert  space  over  the  real  field  R and  E is  a "so-called"  resolution 
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of  identity  (or  spectral  measure)  in  H.  [1]  A resolution  of  identity  E 
is  a family  of  bounded  linear  operators,  E(A),  on  H defined  for  each 
Borel  subset  A of  the  real  number  set  R,  satisfying  the  following  con- 
ditions: 

i)  E(A)  is  an  orthonormal  projection  for  each  A c 8(R)  * 

the  set  of  all  Borel  subsets  of  real  number  set  R;  i.e. 

[EU)]2  » E(A)  * EU)* 

ii)  E(A1)  E(a2)  » E o A2),  V Alt  A2  c 8(R) 

iii)  E*  E((-«»,  t))  is  continuous  in  the  strong  operator  topology, 
goes  to  0 as  T goes  to  -•»,  and  goes  to  1^  (identity  mapping  on  H) 
as  T goes  to  •. 

Defined  as  an  operator  from  (H,  E)  to  (H,  E),  T is  said  to  be  causal 
if 

E*XX  « EtX2=>EtTX1  * E'TXj 

A special  class  of  causal  operators  termed  as  left-miniphase  can  be 
defined  by  the  following  condition 

LtX1  * EtX^==/ EtTX1  a E^TXj 

The  significance  of  being  left -miniphase  is  that  the  inverse  operator 
is  also  causal  when  the  invertibility  of  the  operator  is  guaranteed.  For 
a more  detailed  survey  of  the  properties  of  a Hilbert  resolution  space, 
the  reader  is  referred  to  [1,  14,  20]. 

Hilbert  Resolution  Space  Valued  Random  Variables 

When  talking  about  random  variables  in  a Hilbert  space,  the  structure 
of  the  resolution  of  identity  is  redundant.  But  for  the  purpose  of  this 
paper,  "resolution"  is  included  in  the  title  to  avoid  ambiguity.  Let 
X and  n be  random  variables  [4]  taking  values  in  a Hilbert  space  H with 
appropriate  probability  measure.  X is  said  to  have  finite  first  moment 
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if  for  each  h c H 

E |(X.h)|  < - <A) 

and  E(X,h)  is  continuous  in  h. 

Here  "E"  denotes  the  expected  value  of  a scalar  valued  random 
variable  with  respect  to  the  probability  space  underlying  X.  A random 
variable  satisfying  condition  (A)  has  a "mean"  Mx  which  is  defined  by 
the  equation: 


E(X,h)  « (h,  Mx). 

The  existance  and  the  uniqueness  of  Mx  is  guaranteed  by  condition 
(A)  following  from  the  Riesz  representation  theorem  [21]. 

In  the  sequel  we  deal  with  zero-mean  random  variables  unless  other- 
wise specified.  Such  a random  variable  is  said  to  have  finite  second 
moment  if 

E(X,h)2  < - for  all  h e H (B) 

and  it  is  continuous  in  h. 

As  such,  X has  a "variance  operator"  Qx  which  is  defined  as  follows 
E(X,h)CX,g)  a (h,  Qxg)  V h g c H 

The  existence  and  uniqueness  of  Q is  assured  by  the  representation 
theorem  for  bilinear  functionals  on  Hilbert  space  [9]  together  with 
condition  (B).  More  generally,  given  two  zero-mean  processes  X and  n 
with  finite  second  moment,  a covariance  operator  can  be  defined  as 
follows 

E(X,h)(n,g)  = (h,Q  g)  , V h c H 

AH 

clearly  Qx  a Q^. 

The  covariance  operators  satisfy  the  following: 

U n)  * L QXn  "*•  for  bounded  linear  operators  L and  M 

111  °x  ♦ „ * * Qx„  * 

iii)  QXn  « , in  particular,  Q x * Q* 


iv)  Qx  is  positive,  i.e.  (h,  Q^h)  > 0 
V h c H 

*y 

v)  E| |x| I*  < • if  and  only  if  Q is  nuclear  in  which  case 

A 

E||X||2  * TrCQxJ 

Finally,  we  say  X and  n are  independent  if  * 0. 

II.  The  Optimal  Filter 

With  the  terminology  defined  above,  we  can  now  ftrmally  state  the 
problem  as  follows: 

Let  X and  n be  Hilbert  Resolution  space  (N,E)  valued 
random  ^variables.  X and  n are  zero-mean  and  independent.  X and  n 
have  Qx  and  Qr  as  their  variance  operator,  respectively.  We  want 
to  find  a causal  filter  operator  T on  (H,E)  such  that  the  output 
y of  the  filter  (with  X + n as  the  input)  is  the  best  copy  of  X 
in  the  sense  that  the  error  e (defined  as  X - y)  has  a minimal 
variance  operator  Q#.  Note  here  that  is  a positive  operator  hence 
it  makes  sense  to  talk  about  the  minimal  Q in  the  partial  ordering 
of  positive  operators. 

Even  though  the  problem  is  stated  above  in  its  most  general  form,  it 
is  not  solved  for  the  general  case.  However,  we  do  find  an  optimal  filter 
for  a special  case,  characterized  by  the  following  assumptions: 

i)  Qx  is  Hilbert -Schmidt, 

ii)  Qx  + is  invertible. 

iii)  The  optimal  filter  T to  be  found  is  restricted  to  the  class  of 

Hilbert-Schmidt  operators.  Note:  the  assumptions  are  compatible 
with  the  interpretation  of  X as  a random  signal  and  n as  a noise 
process. 

With  these  assumptions  in  mind,  let  us  derive  the  variance  operator 


» 


> 


» 

■ 


* 


» 


in  IIIWMHI  II  i ii Hi 


« d-f  X - y, 

e x X - y * X - T (X  + n)  =(I  - T)  X + Tn 

Qe  3 Q(I-T)  X + ^Tn 

* ( I-T ) Qx  (I-T)*  + T Q^T* 

3 T (QX  + Qn)  T*  ~ T(3X  " ^ * QX‘ 

In  deriving  Q , independence  of  X and  n has  been  employed.  We  have 
also  noted  that  Q . = Q + Q is  positive  and  self-adjoint.  Following 

A ' D A Q 

the  operator  factorisation  theorem  [7,  19,  20],  there  exists  a Hilbert 

resolution  space*  (H,  E)  and  a linear  operator  F from  (H,  E)  to  (H,  E) 

such  that  (a)  F is  left -miniphase  and  (b)  Q„  + Q = FF*.  Furthermore, 

x n 

F is  invertible  since  Q„  + Q is  invertible.  Therefore  F 7 is  causal. 

x n 

Then  we  can  rewrite  Q as  follows: 

e 


Since 


Hence 


Q = T F F*  T - TQ  - Qv  T*  + Qv 

= [TF  - Qx(F*)'i](P»  T*  - r~1Qx>  + Qx  - Q/*"1  F_1  Q x 
= CTF  - QX(F*)-1]  [TF  - QX(F*)-1]*  + Qx  - Qx  (F  F*)"1  Qx. 

In  the  above  equation,  Q and  Q (F  F*)  ^ Qv  are  both  positive  and  in- 

XX  X 

dependent  of  T.  Therefore,  finding  the  minimum  of  Qe  is  the  same  as 

-l  -1  * 

finding  the  minimum  of  [TF  - Qx  (F*)  ] [TF  - QX(F*)  ] denoted  as  Q(T) 

from  now  on. 

Q(T)  is  in  quadratic  farm.  The  minimum  Q(T)occurs  when  TF  = QX(F*) 

To  fulfill  this  equation,  we  need  a filter  T which  is  equal  to  Qx< F* ) 1 F 1 
* QX(F  F*)’1  = * Qn^'1"  Unfl ortunate ly , this  operator  is  not  necessari- 

ly causal.  In  order  to  find  an  optimal  causal  filter,  we  decompose 
QX(F*)  1 into  two  terms; 

Qx(F*)~i‘  = A t C,  where  A is  the  strictly  anti-causal  part  of 

QX<F*)"1  and  C is  the  causal  part.  The  existence  and  the  uniqueness  of 

*This  space  has  been  shown  to  be  the  reproducing  kernel  resolution  space  for 
Qyt  Q [7,19],  which  does  not,  in  general,  coincide  wi{h.(H,E).  Indeed,  the  factor- 
ization may  not  exist  if  one  requires  that  (H,E)  and  (H,E)  coincide. 

j 


t 
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the  decomposition  are  guaranteed  by  the  fact  that  Qy(F*)~i  is  Hilbert- 
Schmidt.  Readers  are  referred  to  [1,  7,  9,  10,  19,  20]  for  the  termin- 
ology and  details. 

Substituting  QV(F*)”^  back  into  the  equation  for  Q(T),  we  obtain 

A 

Q(T)  = (T  F - C - A)  (T  F - C - A)* 

» (T  F - C)  (T  F - C)*  - (T  F - C)  A*  - A (T  F - C)  + A A* 


In  this  equation,  A(T  F - C)ft  is  strictly  anti-causal  and  (T  F - C ) Aft 
is  strictly  causal. 

By  the  assumption  that  T and  are  Hilbert-Schmidt,  each  term  in  the 
equation  far  Q(T)  is  found  to  be  nuclear.  Therefore,  the  trace  of  Q(T)  is 
taken. 


TrCQ(T)]  = Tr  [(T  F - C)  (T  F - C)*]  t Tr  [A  A*] 

Two  terras  are  dropped  in  the  above  equation  due  to  the  fact  that  the  trace  of 
any  strictly  causal  (or  anti-causal)  nuclear  operator  is  zero  [9,  10].  It 
is  also  known  that  the  trace  of  any  positive  nuclear  operator  is  positive. 

Hence,  the  minimum  of  Tr[Q(T)]  occurs  when  T F = C,  i.e.  when  we  have  a 
filter  T = C F \ This  filter  C F ^ is  causal  following  our  derivation. 

Now  the  only  problem  left  is  to  verify  that  Q(C  F ) is  minimal.  The 
verification  is  straigtforward  following  from  the  fact  that  a positive, 
self-adjoint  operator  is  zero  if  and  only  if  its  trace  is  zero  [20]. 

To  summarize  what  we  have  done  in  this  section,  we  formulate  the 
following  theorem. 

Thm.  i)  Let  X and  n be  Hilbert  Resolution  Space  (H,  E)  valued 
random  variables. 

ii)  X and  n are  zero-mean  and  independent. 

Qv  and  Q denote  the  variance  operators  of  X and  n,  respectively, 
x n 

iii)  is  assumed  to  be  Hilbert-Schmidt  and  is 

invertible. 

iv)  F is  a left-miniphase  factorization  of  Q + Q . 

A II 
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The,  the  optimal  causal  filter  operator  among  the  class  of  the 
HiUbert-Schmidt  operators  is  C F \ where  C is  the  causal  part  of  the 
y operator 

III.  Conclusions 

^ The  Wiener-Hopf  filter  derived  above  is  sketched  diagramatically  in 

Figure  2.  The  first  transformation,  denoted  by  F~\  is  usually  termed 
the  whitening  filter  since  it3  output,  the  innovations  process,  is  white 
I noise;  i.e.  it  has  a memory less  covariance  operator  [7].  It  is  signifi- 

cant that  for  the  filter  to  be  well  defined  it  was  necessary  to  take 

(H,E)  to  be  the  reproducing  kernel  resolution  space  for  Qv  + Q rather 

a □ 

9 than  the  given  resolution  space  (H,E).  As  such,  one  may  conclude  that 

the  innovations  process  naturally  "lives"  in  this  reproducing  kernel 
space  thus  yielding  a further  justification  for  the  study  of  this  ab- 
f stract  resolution  space  even  though  the  given  system  is  defined  on  a 

concrete  resolution  space. 

• 


t 


» 


» 


(H,E) 


Figure  2.  Diagramatic  Representation  of  the 
Wiener-Hopf  Filter 
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Although  the  above  derivation  is  restricted  to  a very  specialized 
filtering  problem  we  believe  that  the  techniques  employed  are  indicative 
of  those  which  will  eventually  lead  to  a general  theory  of  quadratic 
optimization  in  resolution  space.  Indeed,  the  authors  have  already  made 
significant  progress  towards  the  generalization  of  the  above  concepts  to 
stochastic  control  and  estimation  problems  for  systems  described  by  both 
input-output  and  state  models  C20,24].  Of  course,  these  results,  as  with 
the  above  Wiener-Hopf  filter,  apply  to  distributed  and  time-variable 
systems  as  well  as  the  classical  LLF  systems. 

A careful  inspection  of  the  -derivation  of  the  above  filter  will 
reveal  that  the  restriction  that  Qy  + Q be  invertible  can  be  dropped  by 
exploiting  the  fact  that  the  miniphase  factor  of  is  one-to-one 

for  arbitrary  covariance  operators  [7,  19].  As  such,  if  one  replaces  F * 
by  F L (left  inverse)  the  above  derivation  may  be  carried  out  without  the 
invertibility  assumption.  Finally,  we  note  that  the  Hilbert -Schmidt 


assumptions  are  required  only  to  make  the  trace  well  defined.  The  de- 
rivation,  however,  remains  (formally)  valid  if  one  allows  arbitrary  oper- 
ators and  hence  infinite  values  for  the  trace.  Alternatively,  one  may 
replace  the  trace  by  the  memoryless  part  transformator  [1,  10]  obtaining 
the  Wiener-Hopf  filter  by  minimizing  the  memoryless  part  of  the  error 

covariance  (in  the  partial  ordering  of  positive  operators).  This  results 

-R 

is  an  optimal  (though  non-unique)  filter  whenever  the  operator  CQV(F*)  ] 
has  a well  defined  additive  decomposition  into  causal  and  strictly  anti- 

■ 

causal  parts. 

* i : 

: 

j. 
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INTRODUCTION 

Wiener-Hopf  filtering  is  a widely  used  technique  in  certain  kinds  of 
optimization  problems.  The  purpose  of  this  paper  is  to  formulate  Wiener- 
Hopf  filtering  in  abstract  spaces  (reflexive  Banach  resolution  spaces)  and 
to  examine  problems  involved  for  the  formulation  and  the  solving  of  the 
Wiener-Hopf  filter. 

Referring  to  what  has  been  done  in  the  frequency  domain  of  the  classical 
Wiener-Hopf  filtering*,  we've  found  five  major  problems  for  the  formulation 
of  Wiener-Hopf  filtering  in  abstract  spaces.  They  are 

i.  Random  variables  in  abstrace  spaces 

ii.  Causality 

iii.  Operator  factorization 

iv.  Operator  decomposition 

v.  Optimization. 

These  problems  are  briefly  introduced  as  follows: 
i.  Random  process  can  be  thought  of  as  a random  variable  which  takes 
values  in  a function  space.  In  order  to  do  so,  we  need  an  adequate  prob- 
ability measure  over  the  space  involved.  Fortunately,  this  kind  of  prob- 

2 

ability  measure  has  been  defined  over  metric  space  . For  our  purposes, 
we  assume  that  the  space  Involved  is  reflexive  Banach  space,  not  only 
because  this  kind  of  space  possesses  nice  properties  but  also  because 
stochastic  concepts  such  as  "mean"  and  "variance  operation"  can  be  defined 
therein.  Random  variables  taking  values  in  reflexive  Banach  space  is 
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discussed  in  section  II  with  probability  measure  assumed  implicitly. 

f ii.  Concepts  of  causality  have  been  introduced  into  Hilbert  space- 

3 

the  so-called  Hilbert  resolution  space  . In  section  III,  we  extend  the 
works  done  for  Hilbert  space  to  Banach  space.  Concepts  of  causality, 

9 such  as  causal,  anti-causal,  miniphase  and  maxiphase,  are  defined. 

Emphases  are  given  to  reflexive  Banach  resolution  space. 

iii.  Operators  to  be  factorized  in  the  form  of  KK*,  where  K*  de- 

t notes  the  adjoint  of  K,  have  to  be  "positive"  and  "self-adjoint".  These 

commonly- used  properties  among  operators  on  Hilbert  space  can  be  extended 
to  operators  which  map  reflexive  Banach  space  to  its  dual  space.  Factor- 
t ization  theorem  is  given  in  section  IV.  Factor  operator  K is  required 

to  be  left-miniphase. 

iv.  The  decomposition  of  operators  over  Hilbert  space  is  treaded 

* in  Ref.  3.  For  operators  over  Banach  spaces,  this  problem  is  still  under 
research.  For  our  convenience,  operators  are  restricted  to  those  which 
guarantee  the  decomposition. 

♦ v.  As  in  the  classical  Wiener-Hopf  filtering,  we  would  like  to 
minimize  the  variance  of  the  error.  However,  when  Wiener-Hopf  filtering 
is  formulated  in  reflexive  Banach  space,  the  variance  of  the  error  is 

a positive  and  self-adjoint  operator  which  can  only  be  minimized  in  the 
partial  ordering  of  the  positive  operators.  This  subject  is  treaded  in 

section  V. 

I 

BANACH  SPACE  VALUED  RANDOM  VARIABLES 

The  theory  of  Banach  space  valued  random  variables  has  been  studied 
in  Ref.  2.  For  our  purpose,  we  discuss  reflexive  Banach  space  valued 

random  variables  with  probability  measure  over  the  space  assumed  implicitly. 

C 
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The  development  follows  that  of  Parthasarathy  (2)  and  Balakrishnan  (4); 
the  reader  Is  referred  to  these  works  for  the  details. 

Let  p,  ir  denote  finitely  additive  random  variables  taking  values  in 
a reflexive  Banach  space  B.  For  such  random  variables,  we  assume 
E{ | (p , x*) | } < » , for  all  x*  e B* 

(2.1) 

E{(p,  x*J}  is  continuous  In  x* 

Here  E{ * } denotes  the  expected  value  of  a scalar  valued  random  variable 
with  respect  to  the  probability  space  underlying  p.  For  random  variables  which 
satisfy  condition  (2.1),  there  Is  a unique  vector  mp  in  B satisfying 
E{(p,  x*)}  * (nfc,  x*),  for  x*  e B*. 

it*)  is  termed  as  the  mean  of  random  variable  p.  As  in  most  stochastic 
processes,  mean  is  not  our  prime  concern.  Therefore,  in  the  sequel  we 
only  deal  with  zero-mean  random  variables.  For  such  random  variables, 
we  further  assume 


E{  | (p  ,x*)  (tt,  y*)|>  < », 
for  all  x*,  y*  e B* 

(2.2) 

E{(p,  x*)  (*,  y*)> 
is  continuous  in  x*  and  6* 

It  can  be  shown  that  condition  (2.2)  implies  condition  (2.1).  Now 
let's  take  a look  at  E{(p,  x*)  (n,  y*)}.  If  we  fix  y*,  then  E { (p , x*) 

( it , y*)}  is  a bounded  linear  functional  on  B*  (so  an  element  of  B**=B). 
This  means  that  there  exists  a unique  py*  in  B such  that  E{(p,  x*  (it,  y*)} 
a(Py*.  x*)  for  x*  e B*.  Define  a mapping  Qpit  . g*  ^ B>by  Q^*  ,p 

Hence  E{(p,  x*)  (it,  y*)}  ■ (Qpw  y*.  x*) . 
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It  can  be  easily  proven  that  Q Is  linear.  Moreover,  Q Is  bounded. 

P IT  p IT 

t Q Is  termed  as  the  covariance  operator  of  random  variables  p and  n. 

Covariance  operators  satisfy  following  conditions: 

I.  Q(Lp)  ■ L Q M*  , where  L and  M are  linear  bounded  operators 

« on  B. 

II.  Define  Q ■ Q , then  Q ■ Q + Q + Q + Q . Q Is  called  the 

p PP  p^TT  p P7T  Tip  TT  P 

variance  operator  of  p. 

* “ ^np  , In  particular  Qp  - Qp* 

2 

1v.  Qp  Is  positive  In  the  sense  that  (Qp  y*,  y*)  ■ E{(p,  y*)  }_>  0, 
for  all  y*  c B. 

These  conditions  result  from  straight  forward  manipulation  of  the 
defining  equation  for  the  covariance  operator.  Using  Qpii,  we  say  that 
p and  it  are  Independent  If  Q «0. 


BANACH  RESOLUTION  SPACE 

By  a Banach  resolution  space,  we  mean  a 2-tuple,  (B,  gF),  where  B 
is  a Banach  space  and  gF  Is  the  so-called  resolution  of  identity  In  B, 
which  Is  defined  In  the  following: 


Resolution  of  Identit 


Definition  3.1.  Let  B be  a Banach  space.  By  a resolution  of  identity, 
gF,  In  B,  we  mean  a family  of  linear  bounded  operators,  gF  (a),  on  B 
defined  for  each  Borel  subset,  a,  of  the  real  number  set  R,  satisfying 
the  followings: 

1.  gF( R)  ■ Ig-identlty  operator  on  B 
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11.  bfUi)  ' bF(a2)  * bfUi(U2),  for  all  , a2  e 6( R)  - the 

set  of  all  Borel  subsets  of  R. 

n n n 

iii.  dF( Ua  i)  3 I nFUi),  where  {a,}  is  a finite  set  of  disjoint 

B i B 1 

Borel  subsets  of  R. 

iv.  | 1 gF( a ) x 111  ||x||,  for  all  A e 0 (R)  and  x e B (Equivalent 
statement:  Norm  of  gF  (a)  is  either  0 or  1) 

The  subscript  on  the  left  in  the  notation,  gF,  is  to  notify  that 
the  resolution  of  identity  is  defined  over  space  B and  will  be  dropped 
if  no  ambiguity  would  result. 

Working  with  a Banach  resolution  space,  (B,gF),  it  is  natural  and 
important  to  ask  whether  we  can  define  a resolution  of  identity  in  B*, 
the  dual  space  of  B.  The  following  theorem  gives  us  the  answer. 

Theorem  3.1.  Let  (B,  gF)  be  a Banach  resolution  space,  then 
{gF  *(a)  |a  e 0 (a)}  is  a resolution  of  identity  in  B*-termed  as  the 
induced  resolution  space,  (B*,  gF*). 

With  the  resolution  of  identity  defined  as  above,  we'd  like  to  point 
out  that  although  Hilbert  space  is  a special  case  of  Banach  space,  Def- 

3 

inition  3.1  does  not  lead  to  a Hilbert  resolution  space  . In  Hilbert 
resolution  space,  the  resolution  of  identity,  (E(a)  |a  e 0 (R)},  sat- 
isfies an  additional  condition,  i.e.  E*  (a)  = E(a). 

Example.  Let  p,  q e R,  such  that  1/p  + 1/q  3 1.  Then  L is  a reflexive 

r 

Banach  space  with  dual  space  L^.  For  each  f e Lp,  define  (f,q) 

= -f(t)  g(t)  dt.  Let  F (a)  f(t)  3 x( a)  f(t)  =<°  , t/A 

]f( t) ,teA 

It  is  easy  to  show  that  {F  (a)  | Ae0(R)  is  a resolution  of  identity 

in  Lp  and  F*  (a)  * xU)  for  all  a e 0 (R). 
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9 

(B)  Concepts  of  causality 

i Definition  3.2.  Let  (X,XF),  (Y,yF)  be  Banach  resolution  spaces. 

T:  X-*-Y,  is  a linear  bounded  operator 

(i)  T is  causal  if  XF*  ■ XF*  x2  -*•  yF*  T x^  » yF*  T x^,  where 
gF*  * BF(-,t)  , B-  X,  y. 

(ii)  T is  anti-causal,  if  XFt  X1  - XFt  x2  - TFt  T Xj  » yF1  T x2 
where  gFt  = gF  (t,“),  B ■ X,  Y. 

t (iii)  T is  memoryless,  if  T is  causal  and  anti-causal. 

(iv)  T is  left-miniphase  if  xFt  Xj  * xFt  x2  <->  yFt  T x1  3 yFt  T x2 

(v)  T is  left  maxiphase  if  xFt  Xj  * XF^  x2  ^ yFt  T Xj  = yF^  T x2 

C 

(vi)  T is  right-mini  phase,  if  T F^  {X}  = yF^  { Y} 

(vii)  T is  right-maxi  phase,  if  T^FMX}  = yFt  { Y} 

According,  to  above  definitions,  we've  found  the  following  results: 

I 

(1)  Miniphase,  left-or-right-,  implies  causality. 

(2)  Maxiphase,  left-or-right-,  implies  anti -causality. 

(3)  When  X and  Y are  reflexive,  we  have 

9 

(a)  T is  causal  <->  T*  is  anti  causal 

(b)  T is  left-miniphase  <— > T*  is  right-maxi  phase. 

(4)  When  X and  Y are  reflexive  and  T is  invertable,  we  have 
9 

(a)  Miniphases  are  equivalent,  so  are  maxiphases. 

(b)  T is  miniphase  ->  T and  T-1  causal 

T is  maxiphase  — >T  and  T”1  anti-causal. 

Readers  are  referred  to  Ref.  5 for  the  details. 

9 

OPERATOR  FACTORIZATION 

Not  every  operator  over  arbitrary  Banach  spaces  can  be  factorized 
in  the  form  desired.  For  our  purpose  the  desired  form  of  factorization  is  K K*. 
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An  operator  to  be  factorized  in  this  form  has  to  be  positive  and  self-adjoint. 
These  two  commonly-used  properties  for  operators  over  Hilbert  space  can  be 
extended  to  operators  that  map  from  reflexive  Banach  space  to  its  dual 
space.  They  are  defined  as  follows: 

Definition  4.1. 

(i)  Let  B be  a reflexive  Banach  space. 

(ii)  Q=  B-»B*,  is  linear  and  bounded.  Q is  said  to  be  positive  if 
(x,  Qx)  _>  0,  for  each  x e B.  Q is  said  to  be  self-adjoint  if  Q*  s Q. 

Note  that  Q*  : B**=B  -*■  B*  so  it  makes  sense  to  compare  Q with  Q*. 

For  positive  and  self-adjoint  operators,  we  have  the  following  theorem: 
Theorem  4.1. 

(i)  Let  B be  a reflexive  Banach  space. 

(ii)  Q : B-*8*,  is  linear,  bounded,  positive  and  self-adjoint. 

Then  there  exist  a Hilbert  space  H and  a linear  bounded  operator  K : 

H-*B , such  that  Q s K K*®. 

When  dealing  with  Banach  resolution  spaces,  the  usefulness  of  operator 
factorization  is  limited  unless  the  factor  operator  possesses  certain  causal 
properties.  Referring  to  factorization  of  the  spectral  density  in  classical 
Wiener-Hopf  filtering,  we  have  found  what  we  need  is  a factorization  theorem 
which  gives  a causal  operator  and  guarantees  a causal  inverse  once  the 
existance  of  a inverse  is  granted,  i.e.  a theorem  that  gives  a miniphase 
factorization.  Based  on  Theorem  4.1,  we  construct  the  resolution  of 
identities  in  spaces  involved  and  we  come  up  with  the  following  theorem. 
Theorem  4.2. 

(i)  Let  (B,F)  be  a reflexive  Banach  resolution  space.  (B*.  F*) 
denotes  the  induced  resolution  space. 
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(ii)  Q * (B,F)  -►  (B*.  F*) , is  linear,  bounded,  positive  and  self- 
adjoint.  Then  there  exist  a Hilbert  resolution  space  (H,E)  and  a linear 
bounded  operator  K * (H,  E)  -*■  (B*,  F*)  such  that 

1.  Q - K K* 

2.  K is  a left-mini  phase 

3.  The  factorization  is  unique  up  to  a memoryless  unitary 
transformation. 

For  the  proof  of  this  theory,  please  refer  to  Ref.  5. 

HIENER-HOPF  FILTERING  FORMULATED  IN  REFLEXIVE  BANACH  SPACE 

With  the  preparation  of  sections  II,  III  and  IV,  now  we  are  ready 
for  the  formulation  of  Wiener-Hopf  filtering.  The  formulation  is  done 
as  follows: 

Let  X,  n be  random  variables  taking  values  in  a reflexive  Banach 
resolution  space  (B,  F).  X denotes  the  signal  and  n the  noise.  Both 
X and  n satisfy  condition  (2.2)  in  section  II  and  they  are  assumed  to 
be  zero-mean  and  independent.  As  such,  X and  n have  Qx  and  Qn  as  their 
variance  operators  respectively.  The  problem  we  are  facing  is  to  find 
a filter,  T : B-*8,  linear  and  causal,  to  operate  on  X + n such  that  the 

error,  defined  as  x-y,  where  y is  the  output  of  T,  has  a variance  operator 
that  is  minimal  in  the  partial  ordering  of  the  positive  operators.  We 
will  describe  this  ordering  right  after  we  find  the  variance  of  the  error. 

Let  e denote  the  error  and  Q denote  its  variance  operator.  Since 
def  e 

e a X-y  * X - T(X+n)  ■ (I-T)X  +Tn  we  have 

Qe  a ( I -T)  Qx  ( I-T)*  + T Qn  T*,  following  from  the  results  in  section 
II.  Rearranging  terms  in  Q , we  get 
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Qe  = T(QX+Qn)T*  " Qx  T*  " T QX  + QX‘ 

Qg  is  dependent  on  T.  We  write  Qg  (T)  to  notify  the  dependence. 

CL  (O  is  said  to  be  minimal  for  some  filter  T , if 
e o o 

Qe(T)  £Qe(T0)  ->  Qe(T)  « Qe(T0)  (A  < B if  (B-A)  is  positive). 

In  the  equation  for  Q , Q + Q represents  the  variance  operator 

of  X+n,  hence  is  positive  and  self-adjoint.  Therefore,  by  Theorem  4.2, 

there  exists  a resolution  space  (H,E)  and  a linear  bounded  operator 

K = (H,E)  -*■  (B,F) , such  that  (a)  Qg  = K K*  (b)  K is  a left-miniphase. 

Without  further  assumptions,  the  formulation  would  be  stuck  right  here. 

At  this  point,  what  we  need  is  an  invertable  factor  operator  K.  The 

invertability  of  K can  be  guaranteed  by  the  invertability  of  Qx  + Qn> 

There  are  several  ways  to  secure  the  invertability  of  Qx  + Qp.  One  way 

is  to  assume  that  Qx  + Qn  in  onto  and  Qn  is  positive  definite.  With  an 

invertable  factor  K,  can  be  rewritten  an 

e 

Qe  = T K K*  T*  - Qx  T*  - K Qx  + Qx  = { TK-QX( K*) " 1 } {K*T*-T-1QX>  + 

Qx  - QX(K*)_1  K'1  Qx  = (TK-QX(K*)'1}  { TK-QX( K*) " 1 } * + Qx  - Qx  (KK*)_1QX 

The  last  two  terms  in  the  above  equation,  Qx  and  QX(KK*)~^QX,  are 

positive  and  independent  of  T.  Hence,  to  find  the  minimal  of  Qg  is  the 
same  as  to  find  the  minimal  of  {TK-QX(K*)”*}  { TK-Q^  (K*)”*}*  - denoted  as 
Q(T)  in  the  sequel.  Right  now,  we  are  facing  the  same  kind  of  problem  as 
in  classical  Wiener-Hopf  filtering.  Minimal  Q(t)  occurs  when  T = QX(K*)~*K~*, 
but  it  does  not  represent  a causal  system  in  general.  In  order  to  get  a 
possible  optimal  causal  filter,  can  we  decompose  QX(K*)-1  into  "causal  part" 
and  "strictly  anti-causal"  (a  term  to  be  generalized  in  Banach  resolution 
space)  and  under  what  conditions  can  we  do  so?  This  subject  has  been 
treated  in  Ref.  3 and  Ref.  7 when  the  reflexive  Banach  space  happens  to 
be  a Hilbert  space.  However  in  reflexive  Banach  resolution  spaces,  the 
subject  is  still  under  research.  While  we  follow  the  same  pattern  as  that 
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of  classical  Wlener-Hopf  filtering  in  frequency  domain,  we  would  like 
to  ask  whether  this  decomposition  would  work  and  how  It  would.  The 
same  question  in  classical  Wlener-Hopf  filtering  Is  not  directly  answered 
In  frequency  domain.  In  order  to  find  the  answer,  let's  assume  the 
decomposition.  Let 

Qx(K-)-1  ■ C + A,  where  C Is  the  causal  part  of  QX(K*)”*  and  A Is 
the  "strictly  anti-causal  part"  (a  term  to  be  generalized  In  Banach 
resolution  space).  Then 

Q(T)  - (TK-C-A1  {TK-C-AJ* 

■ { TK-C > iTK-Cl*  - A{ TK-C}* 

» (TK-C)A*  + AA* 

To  claim  TK-OO  Is  the  condition  for  minimal  A(T),  we  should  demon- 
strate that  those  cross  terms, (TK-C)  A*  and  A{TK-C}*,  have  no  effect  on 
the  ordering  of  Q(T).  Again  when  the  reflexive  Banach  resolution  space 
Is  Hilbert  resolution  space,  we've  found  two  ways  to  achieve  this.  The 
first  one  is  to  take  the  trace  of  Q(T).  Surely,  work  has  to  be  done  to 
guarantee  Q(T)  being  nuclear.  The  second  one  is  to  take  the  memoryless 
part  of  Q(T).  This  is  justified  once  the  decomposition  is  given.  However 
we've  also  found  advantages  and  disadvantages  to  each  way.  For  the  method 
of  taking  trace,  it  gives  a minimal  variance  operator  once  the  maximum 
of  the  trace  is  found,  but  we  have  to  restrict  certain  operators,  such 
as  Qx  and  T,  to  be  Hi lbert-Schmidt.  On  the  other  hand,  the  method  of 
taking  memoryless  part  works  for  a broader  class  of  operators-operators  which 
have  decompositions,  but  it  does  not  give  a minimal  variance  operator. 

The  best  we  can  have  is  a variance  operator  that  has  a minimal  memoryless 
part.  However,  there  is  an  important  aspect  for  taking  the  memoryless 
part.  This  method  allows  us  to  generalize  the  idea  in  reflexive  3anach 
space,  while  the  other  method  does  not.  The  reason  is  quite  simple,  for 
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It  does  not  make  sense  to  talk  about  the  eigenvalue  of  an  operator  that 
maps  from  Banach  space  to  Its  dual,  not  to  mention  the  trace  of  such  an 
operator,  while  It  does  not  make  sense  to  take  the  memoryless  part  given 
the  decomposition.  Readers  are  referred  to  Ref.  5 for  the  details  of 
Wlener-Hopf  filtering  formulated  In  Hilbert  resolution  space.  When  all 
the  problems  mentioned  above  are  solved,  we  would  come  up  with  the  optimal 
filter  T0  ■ C K"1,  a causal  system. 


CONCLUSION 

Wiener-Hopf  filtering  has  been  formulated  and  solved  In  Hilbert 
resolution  space5.  In  this  paper,  we  outlined  the  formulation  In  re- 
flexive Banach  resolution  space  and  the  possible  way  of  solving  It. 
Generalization  would  be  accomplished  once  the  theory  of  operator  de- 
composition in  3anach  resolution  space  is  completed. 
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Reproducing  Kernel  Resolution  Space  and  Its 
Applications  II* 
l.  Tung*  and  R.  Saeks0 

Abstract 

This  paper  extends  the  concept  of  a reproducing  kernel  resolution  space  to  a 
Banach  space  setting.  The  resultant  reproducing  kernel  resolution  space,  however, 
remains  a Hilbert  space  thereby  permitting  a number  of  problems  In  mathematical 
system  theory  to  be  transformed  from  Banach  space  to  Hilbert  space.  Particular 
emphasis  Is  given  to  the  study  of  Banach  space  valued  random  variables  and  the 
scattering  operator  formalism  for  an  electric  network. 


I.  Introduction 

12 

In  a previous  paper  one  of  the  authors  exhibited  the  relatlonsnlp  between 
the  factorization  problems  which  arise  In  mathematical  system  theory  and  the 
reproducing  kernel  resolution  spaces  introduced  by  Kallath  and  Outtweiller^.  The 
purpose  of  the  present  paper  Is  to  show  that  much  of  this  work  can  be  extended  to 
a Banach  space  setting  without  the  loss  of  Its  Hilbert  space  character.  Indeed, 

It  Is  shown  that  the  reproducing  kernel  resolution  space  for  a positive  self-ad- 
joint  operator  mapping  a reflexive  Banach  space  to  its  dual  is  a Hilbert  space. 
Since  this  Is  precisely  the  class  of  operators  encountered  when  the  factorization 
problems  of  mathematical  system  theory  are  extended  to  a 3anach  space  setting,  the 
resultant  theory  allows  one  to  transform  systems  problems  from  a 3anach  space  to 
a Hilbert  space  setting.  In  particular.  It  Is  shown  that  the  study  of  certain 
Banach  space  valued  random  variables  can  be  carried  out  in  the  reproducing  kemal 
Hilbert  resolution  space  defined  by  Its  covariance.  Secondly,  it  is  shown  that  an 
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electric  network  with  voltage  and  current  vectors  defined  In  Banach  space  may  be 
characterized  by  a scattering  operator  defined  on  an  appropriate  reproducing 
kernel  Hilbert  resolution  space. 

A Banach  resolution  space  Is  defined  and  Its  elementary  properties  developed. 

We  note  that  the  axioms  for  a Banach  resolution  space  are  weaker  than  those  re- 
quired fora  Hilbert  resolution  space  and,  as  such,  the  theory  developed  does  not 
necessarily  specialize  to  the  classical  Hilbert  resolution  space  theory.  The 
axioms  are,  however,  sufficient  for  the  present  purposes.  In  particular,  such  con- 
cepts as  causal,  antlcausal,  miniphase,  and  maxiphase  operators  are  well  defined. 

In  the  third  section  a factorization  theorem  for  operators  mapping  a reflex- 

6 7 3 

Ive  Banach  space  to  Its  dual  recently  developed  by  Chobanlan  , Vakhanla  , and 

g 

Masanl  Is  applied  to  Banach  resolution  spaces  to  develop  miniphase  and  maxiphase 
factorization  theorems.  With  the  help  of  these  theorems  a "unique"  factor  space 
- the  generalization  of  the  RKRS  to  Banach  resolution  space  - Is  formulated. 

While  sections  II  and  III  deal  with  fundamental  theorems,  sections  IV  and  7 
are  devoted  to  the  application  of  these  theorems.  The  first  application  con- 
sidered Is  the  study  of  reflexive  Banach  space  valued  random  variables.  With  the 
probability  measure  over  the  Banach  space  assumed  Implicitly  stochastic  concepts 
such  as  mean  and  covariance  operator  are  defined.  Since  the  resultant  variance 
operator  Is  a positive  self-adjoint  mapping  from  a reflexive  Banach  space  to  its 
dual  the  factorization  theory  developed  In  the  previous  section  may  be  Invoked 
to  transform  the  given  random  variable  Into  Its  reproducing  kernel  Hilbert  resolut- 
ion space. 

In  section  V similar  factorization  techniques  are  used  to  define  the  scattering 
variables  for  an  electric  network  characterized  by  voltage  and  current  vectors 
taking  values  In  a reflexive  Banach  space  and  Its  dual,  respectively.  Here,  the 
normalizing  operators  take  the  form  of  maps  from  the  given  3anach  space  to  an 
» appropriate  reproducing  kernel  Hilbert  resolution  space  resulting  In  a scattering 
operator  which  is  defined  on  a Hilbert  space. 
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II.  Banach  Resolution  Space 

By  a Banach  resolution  space,  we  mean  a 2-tuple,  (3,gF),  where  B Is  a Banach 
spaca  and  gF  Is  the  so-called  resolution  of  Identity  In  B,  which  Is  defined  as 
follows. 

Resolution  of  Identity  In  Banach  Space 

Def.  II. 1.  Let  3 be  a Banach  space.  By  a resolution  of  Identity,  gF, 
In  B,  we  mean  a family  of  bounded  linear  operators,  gF(a),  on  3 defined  for 
each  Borel  subset,  a,  of  the  real  number  set  R,  satisfying  the  following 
condl tlons: 

1.  gF(R)  ■ Ig  ■ Identity  operator  on  B. 

11.  gPCa-j)  8F(a2)  • gF(a1  (1  a2),  for  all  ^ , a2  c a(R)  ■ the 

set  of  all  Borel  subsets  of  R. 
n n n 

111.  gF(U  A.j)  ■ c gF(a^) ,{a^ : finite  set  of  disjoint  Borel 
subsets  of  R. 


1v.  ||BF(a)x||  ^ |I*||,  for  all  a c s(R)  and  x c B (equivalently 

the  norm  ofgF(a)  Is  either  0 or  1). 


The  subscript  on  the  left  In  the  notation,  gF,  Is  to  signify  that 
the  resolution  of  Identity  Is  defined  In  space  B,  and  will  be  dropped  If  no  am- 
biguity would  result. 

Working  with  a Banach  resolution  space  (B,  gF) , It  Is  natural  and 
Important  to  ask  whether  we  can  define  a resol uton  of  Identity  In  3*,  the  dual 
space  of  B.  The  following  theorem  gives  us  the  answer. 

Thm.  II. 1.  Let  (B,  gF)  be  a Banach  resolution  space.  Then 
(gF(a)*|a  c a(R) } Is  a resolution  of  Identity  in  B*.  This  resolution  of  Iden- 
tity In  3*  Is  called  the  Induced  resolution  of  Identity. 

The  proof  of  this  theorem  follows  from  straightforward  manipulation  of  the 
definitions  of  adjoint  operator  and  resolution  of  identity. 

The  concept  of  a resolution  of  Identity  can  be  best  understood  by 


r 


31 


examples.  Two  typical  examples  are  Illustrated  as  follows. 

Example  1.  lp:  the  Banach  space  of  equivalence  classes  of  functions 
that  map  R to  R and  satisfy  the  following  Inequality, 

< _ 

£ |f(t)|p  dt  < where  1 < p < •. 

Oeflne 

0 , for  t i a 

FU)f(t)  • xU)f(t)  •<  , for  f c L_. 

jf(t),  for  t « a p 

It  Is  easy  to  show  that  {F(a)|a  c 8(R)>  Is  a resolution  of  identity  In  Lp  once 
the  properties  of  xU)  are  explored. 

Example  2.  Let  p,  q c R such  that  1/p  ♦ 1/q  ■ 1.  Then  Lp  Is  a re- 
flexive Banach  space  whose  dual  space  is  L^.  For  all  f c Lp,  g c L^, 
define 

(f.g)  ■ £ f(t)g(t)  dt. 

Let  {F(a)|a  c a(R) > be  the  resolution  of  identity  in  Lp  as  defined  In  Example  1. 
Then  we  have 

(F(a)f,  9)  * CxU)f(t)]g(t)  dt 
• f(t)g(t)  dt 

" C f(t)[x(a)g(t)]  dt 
■ (f.  xU)g) 

So  {x(a)|a  * 3(R)>  Is  also  the  Induced  resolution  of  Identity  in  Lp. 


m ' 

Y*  ' * 

11.  T Is  said  to  be  anti-causal  If  t 

XFt  X1  * XFt  x2  YFt  T X1  “ YFt  T x2 » w*»«re  BFt  " V BF  , B ■ X,  Y. 

111.  T Is  said  to  be  memory less  If  T Is  causal  and  anti-causal. 

1v.  T Is  said  to  be  left-mini  phase  If 

yF*  T x ■ 0 <■*  yFc  x • 0, 

« 

v.  T Is  said  to  be  left-maxi  phase  If 


yF^  T x * 0 <■*  x ■ 0, 

vl.  T Is  said  to  be  right-mini  phase  If 
T[xFtCX]]  - yFt[Y], 

vl 1 . T Is  said  to  be  right -maxi phase  If 
TC/^X]]  - yFt[Y]. 


i 


Based  on  the  above  definitions,  the  following  results  may  be  readily  verified? 

1.  If  T is  mini  phase,  left-  or  right-,  then  T is  causal. 

2.  If  T Is  .maxi  phase,  left-  or  right-,  then  T Is  anti-causal. 

3.  When  the  spaces  Involved  are  reflexive  and  T X * Y Is  bounded  and  linear,  then 

(a)  T Is  causal  Iff  T*  is  anti-causal 

(the  Induced  resolutions  of  identity  have  been  used 
for  the  determination  of  the  causality  of  T*) 

(b)  T Is  1 eft -mi ni phase  iff  T*  Is  right-maxi  phase 

4.  When  the  spaces  Involved  are  reflexive  and  T:X  * Y-.is  invertible. 

(a)  T Is  left-mini  phase (maxi  phase)  If  it  is  right -mini  phase  (maxiphase). 

(b)  T is  miniphase  ->  T * T“l  are  causal 

T Is  maxiphase  ->  T * T”  are  anti -causal. 

The  reader  Is  referred  to  Ref.  2 for  the  details. 

III.  Operator  Factorization  and  RKRS 

A number  of  applications  arise  in  system  theory  wherein  it  is  desired  to 

factor  an  operator,  Q,  either  In  the  form  K1C*  or  T*T.  For  an  operator  to  be 
factorized  in  either  of  these  forms  it  has  to  be  "positive'*  and  "self-adjoint". 

Although  undefined  In  general,  these  two  commonly  used  properties  for  operators 
over  Hilbert  space  can  be  extended  to  operators  that  map  a reflexive  Banach  space 
to  Its  dual  space.  They  are  defined  as  follows: 


. 
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Oaf.  Ill . 1 . 1.  8 Is  « reflexive  Banach  space. 

It.  T:  B ♦ 8*.  Is  linear  and  bounded. 

T Is  said  to  be  positive  If 

(x,  T x)  » 0,  for  all  x c B. 

T Is  said  to  be  self-adjoint  If 
T*  - T. 

Note  that  T*:  B**»8  8*,  so  that  It  makes  sense  to  equate  T and  T*. 

For  positive  and  self-adjoint  operators,  we  have  the  following  theorem.  'The 
theorem  Is  stated  without  proof.  Interested  readers  are  referred  to  Masanl's  work 
(9). 

Thm.  III.l.  Let  Q:  B -*■  B*»  be  a linear,  bounded,  positive  and  self-adjoint  oper- 
ator, where  B is  a reflexive  Banach  space  and  8*  Is  Its  dual.  Then  there  exists  a 
Hilbert  space  H and  a linear  bounded  operator  X mapping  H to  8*  such  that 
Q • X X*.6, 7,8,9 
Left-  and  Right-  Factorization 

Although  Thm.  III.l  yield  a factorisation  through  a Hilbert  space,  which  Is  valid 
for  any  linear,  bounded,  positive,  self  adjoint  operator  mapping  a reflexive  Banach 
space  to  It  dual  our  applications  In  System  Theory  require  that  the  factors,  K,  have 
appropriate  causality  properties.  Based  on  Thm. III.l,  we  construct  appropriate  reso- 
lutions of  Identity  In  the  spaces  Involved  which  yield  the  following  theorems. 

Thm. I II. 2(1  eft-factorization) . 

1.  (3,  F)  Is  a reflexive  Banach  resolution  space. 

11.  Q:  (B,  F)  -*«  (B*,  F*),  where  F*  Is  the  Induced  resolution  of  identity 

In  3*.  Q Is  linear,  bounded,  positive  and  sel f-adjolnt. 

Then  there  exists  a Hilbert  resolution  space  (H,E}  + and  a linear  bounded 
« operator  K:  (H,  E)  (B»,  F*) , such  that 

1.  Q • K X*, 

2.  X Is  left-mini  phase, 

• 3.  The  factorization  Is  unique  up  to  a memory less  unitary  transformation. 

t The  requirements  for  a Hilbert  resolution  space  are  more  restrictive  than  those  for 
a Banach  resolution  space,  namely  E(a)*  • E(a)  Is  required  for  a Hilbert  resolution 
space,  which  condition  does  not  make  sense  in  a Banach  resolution  space. 
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The. Ill . 3(r1ght- factorization) . 

1.  Q:  (B,  F)  - (B*,  F*). 

11.  Q Is  linear,  bounded,  positive  and  self-adjoint.  Then  there  exists 
a Hilbert  resolution  space  (Jj,  £)  and  a linear  bounded  operator  T: 

(8,  F)  * (JJ,  |),  such  that 

1.  Q • T*  T, 

2.  T Is  right-mini  phase,' 

3.  The  factorization  Is  unique  up  to  a memoryless  unitary  transformation. 

The  proofs  for  the  above  theorems  are  time-consuming  and  they  are  not  our  prime 

concern.  Hence,  they  will  not  be  presented  here.  Interested  readers  are  re- 

ferred to  the  Appendix  for  their  proofs. 

Reproducing  Kernel  Resolution  Space  (RKRS) 

There  is  a conwon  statement  In  each  theorem  of  the  previous  paragraph,  i.e. 

"The  factorization  Is  unique  up  to  a memoryless  unitary  transformation- . For  certain 
applications  such  as  the  study  of  Banach  space  valued  random  variable.  We  would  like 
to  eliminate  this  ambiguity.  This  Is  achieved  via  the  concept  of  a reproducing 
kernel  resolution  space.  First,  wi  define  some  notation. 

Def. III. 2.  Q,  K,  ant  T are  defined  as  in  the  previous  paragraph 

i.e.  Q ■ K !C*  ■ T*  T.  Let  Hg  - R(X)t,  Hg  - R(T*)t.  For  x,  y In  R(K) , let 

A* 

(x,  y)  Hg  ■ (K'L  x,  y)H,  and  for  w,  z In  R(T*),  let 
(z.  w)Hq  » (T*‘L  z.  T*~l  w)  H. 

Define  gE*  ■ K E1  K*L  on  the  Hg  and  Q*t  * T*  E^T*”*"  on  Hg,  then  extend  them  to 
s(R).  Mote  that  E corresponds  to  E in  Thra.III.2. 


R(X)  denotes  the  range  of  the  operator,  X. 


Then,  we  have  the  following  theorem. 

Thn.III.4.  (Hq,  q£)  (the  so-called  RKRS)  defined  above  Is  a Hilbert  resolution 
space  which  Is  Independent  of  the  factorization  Q • K K*  used  In  Its  definition. 
Moreover, 

“V  ^ 

1.  (Hq,  qE)  Is  unltarlly  equivalent  to  (H,  E) , 

11.  R(Q)  Is  dense  In  Hq, 

% % 

111.  (fV  x • 0 Iff  E*  x • 0,  for  x c Hg, 

1v.  K:  (H,  E)  •*  (Hq,  qE),  Is  memoryless, 

where  (H,  E)  corresponds  to  (H,  E)  In  Thm.  III.  2. 

% % 

Proof:  1.  Hq  ■ R(K),  so  Hq  Is  a linear  vector  space. 

11.  By  definition,  (x,  y)jj  ■ (K"L  x,  k"L  y)Jj[.  It  Is 
trivial  to  show  that 

% 

1.  (x,  y)  ■ (y,  x),  for  all  x,  y c Hq,. 

2.  (x,  ay)  • a(x,  y),  for  all  a c R,  x,  y t HR 

3.  (x,  y ♦ z)  • (x,  y) ,+>  (x.  z)  for  all  x,  y,  z c Hg, 

. . “V 

4.  (x,  x)^J  ■ (K"L  x,  1C”1"  x)JJ|  >,  O.for  all  x c Hq  and 

since  K Is  linear,  x /*  0 Implies  K*L  x t 0.  So 

(x,  x) > 0,  for  x c Hq  and  x f 0. 

So  (x,  y)H  an  innep  pro<iUc.t  over  u 

W ^ ^ -L 

111.  Let  (x^>  be  a Cauchy  sequence  In  Hq,  then  (K  x^}  Is 

Cauchy  In  H.  So  x^  * z,  for  some  z c H.  3ut 

. A#  % 

z * P I z,  so  X|  * K : In  H.  Hence  Hg  Is  a Hilbert 
space. 

1v.  Assume  Q • X K+  • K'  K'*,  both  K and  X'  are  left-minipnase 

A,  A,  % \» 

factorizations  of  Qon  factor  spaces  (H,  E)  and  (H * , S'),  re- 
spectively. 

1.  By  Thm. III. 2,  there  exists  a memoryless  unitary  transfor- 
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nation  U:  (H.  if)  - (H\  if'),  such  that  K U • K*.  So 
R(K')  ■ R(K  U)  • R(K) , slnca  U Is  onto.  Hence  ifq  Is 
Independent  of  the  factorization. 

2.  (X’*L  z,  <*’L  w)pj,  • ((X  U)‘L  z,  (X  U)‘L 

- (U*1  X*L  z.  U"1  X-L  w)fl, 

■ (X-**  z,  X'L  w)p|,  for  all 

z,  w c fifq.  So  the  Inner  product  Is  Independent  of  the 
factorization. 

3.  K'  l'1  X’"L  • (X  U)  if’*  (X  U)‘L  -XU  U"1  X'L 

• X Efc  U U']  rL  - X Et  X-1 

Therefore,  gE  Is  also  Independent  of  the  factorization. 

v.  Define  X:  ft  ifq  by  £ x ■ X x,  for  all  x « ft.  Then  £ Is  1-1 

% 

and  onto  and  for  all  x « H, 

II*  xli?  " ||X  x||*  " I Ik’L  < x||?  « 11*11* 

Hg  rfg  ft  H 

So  X Is  a unitary  mapping.  Furthermore,  we  have 

gE  • X E‘  X on  Hg  and  for  all  z c Hg,  there  is  a unique 

x c if  such  that  z • X x - X x.  So  x « X‘L  z - X*1  z, 
f.e.  X u • X . Hence  gE1*  ■ X E*  X . This  means  that 
(Hg,  gE)  is  unltarlly  equivalent  to  (H,  E). 

vl . R(Q)  ■ R(X  X*)  - X[R(K*)]  • ft[R(X*)] 

Since  R(X*)  is  dense  in  H,  so  X[R(X*)]  is  dense  In  Hg 
(for  X Is  unitary), 
vll.  1.  X is  1 eft -flrinl phase,  so 


E“  x ■ 0 iff  (FV  X x • 0,  for  x ‘"'H.' 

2.  For  all  z « Hg,  there  exists  x « ft  such  that 
z • X x or  x ■ X’**  z.  So 


1 


J7 


E*  <‘L  Z • 0 Iff  (fV  z • 0,  z * Hg. 

3.  If  (FV  2-0,  then 

: ■ x ? I*L  : « 0. 

And  If  gE*  z • 0,  l.e.  K E*  lx'1*  z • 0,  then 

E*  K"L  z • 0.  So  (FV  z » 0. 
vlll.  K:  (H,  E)  ♦ (Hg,  g£) , Is  actually  K defined  In  v.  K 1$  left-mlnlphasa, 
so  E*  x - 0 Iff  (Ffc)*  K x - 0.  Hence  E*  x • 0 Iff  g£*  K x » 0,  for 
x c H.  But  since  K ■ K,  the  above  equation  Implies  K Is  also  a 
left -mini phase.  So  K*  * K Is  causal  (by  being  left-miniphase) , 

*>.  ■r. 

and  < Is  also  anti -causal.  Therefore,  K Is  memoryless.  i 

Following  Immediately  from  the  previous  theorem  Is  a rather  Interesting  result  which 
gives  us  a sort  of  "unique"  left-factorization.  This  result  is  indicated  In  the  next 
corollary. 

Cor.  Let  (Hg.gE)  be  defined  as  In  Thm.III. 4.  Then  there  exists  a left- factorization 
P:  (Hg,  gE)  * (B*.  F*),  such  that 

1.  P z ■ z for  all  z 1 Hg  c B*, 

11.  P*  b * Q b,  for  all  b « B. 

Proof:  By  Thm.III. 4,  we  have  <:  H - Hq,  a memoryless  unitary  operator,  and 
5!  x » K x,  for  all  x c fif.  Define 

P • K 5T1 : Hg  * B*.  Then 

1.  For  all  z * Hg,  there  exists  x c h such  that  £ x • z. 

2.  By  Thm.III. 4, 

gE1  z ■ 0 Iff  (F1)*  z • 0,  for  z c Hg.  So 

qE*  z ■ 0 Iff. (FC)*  ? z • 0,  for  z c Hq.  Hence  P is 

left -mini  phase. 
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3.  P*  b ■ (K  K’1)*  b ■ (K“V  <*  b • K K*  b ■ K K*  b ■ Q b, 
for  all  b c B.  And  P P*  b ■ P (q  b)  - q b.  # 

Tw  *>/ 

The  ''uniqueness1'  we  mentioned  Is  due  to  the  fact  that  (Hg,  gE)  Is  Independent 
of  the  factorisation. 

As  In  the  previous  paragraph,  there  are  corresponding  dual  theorems  to  Thm.  III. 4. 
These  theorems  are  described  below,  with  proofs  only  sketched. 

Thm. III. 5.  (Hg,  q^) , defined  at  the  beginning  of  this  section.  Is  a Hilbert 
space  which  Is  Independent  of  the  factorization,  q • T*  T,  used  In  Its  definition. 
Moreover , 

1.  (jtfq.  g§)  Is  unltarlly  equivalent  to  (H,  £), 

11.  R(q)  Is  dense  In  Hg, 

111.  (F4)*  x ■ 0 Iff  Efc  x • 0,  for  x c H, 

T’:  <«•  V * («q-  g£) , Is  memoryless.. 

Proof:  1.  The  proof  of  (Hg,  gE)  being  a Hilbert  resolution  space 

and  being  Independent  of  the  factorization  Is  essentially 
the  same  as  that  of  Thm. I I I. 4,  and  therefore  Is  omitted. 

11.  Define  T*  : (H,  E)  (Hq,  g£) , by  T*  x - T*  x for  all 
x c |1.  Then  T Is  1-1  and  onto.  For  all  x c jl, 

lli1*  xll2  • I IT*  x||2  -II  (T*)-*-  T*  x|L2  ■ ||x||J. 

So  T Is  a unitary  mapping.  With  T defined  as  above.  It 
Is  routine  to  verify  the  rest  of  the  theorem.  But  we  need 
to  note  that  T*:  (jl,  £)  - (y~,  g£)  In  1v.,  Is  actually 
the  defined  above. 


Cop.  (JJq,  q£)  defined  as  above,  then  there  Is  a right-factorization  (mini- 
phase) of  Q,  q:  (8,  F)  * ({jq.  g£) , such  that 

1.  q x • Q x,  for  all  x c B, 

11.  q*  z ■ z,  for  all  z * Jlq  c 8*  . 

Proof:  By  Thm.III.5,  |J  * (Jg,  Is  memoryless  and  unitary. 

Define  q - h T:  (B,  F)  - (j^,  <£) . Then 

1.  q x-T*  T x » T*  (T  x)  • Q x,  for  all  x c 8. 

11.  q*  z - (T*  T)*  z • T*  T z ■ T*(T  z)  - z,  for  all 

% 

z c )1q,  since  T*  Is  unitary. 

111.  By  Thm.III.5, 

<&  x " 0 1ff  (Ffc)*  x ■ 0,  for  x c (1.  So 

Q^t  x ■ 0 Iff  (Ffc)*  q*  x ■ 0,  for  x < {1.  So  q*  Is 
left-maxi phase,  l.e.  q Is  rlght-mlnlphase. 

1v.  q*  q b • (?*  T)*  ft*  T)  b - T*  } }*  T b - T*  T b - Q b, 
for  all  b c B. 


One  way  to  view  a random  process  Is  to  consider  It  as  a random  variable  which 
takes  values  In  a function  space.  Of  course,  we  have  to  use  an  adequate  proba- 
bility measure  to  make  the  Idea  work.  Fortunately,  this  kind  of  measure  has 
been  defined  for  metric  spaces.*’*  In  this  section,  we  first  define  stochastic 
properties  such  as  "mean”  and  "variance  operator"  for  a Banach  space  valued 
random  variable.  We  then  look  Into  the  factorization  of  the  variance  operator 
and  the  results  that  can  be  derived  therefrom,  l.e.  the  RKRS.  Interestingly 
enough,  the  RKRS  of  a Banach  space  valued  random  variable  is  a Hilbert  space. 

This  seems  to  be  a nice  result,  but  there  are  obstacles  for  further  application. 
All  these  will  be  discussed  In  the  following. 


(1)  CovarlanceJJgerator 

A probability  measure  on  Banach  space  Is  a rather  complicated  matter.  In  the  sequel, 
we  Implicitly  assume  Its  existence  as  Indicated  by  the  expected  value  symbol  £{•}.  The 
reader  Is  referred  to  reference  4 for  the  details. 

Let  p,  * denote  finitely  additive  random  variables  taking  values  In  a relexlve 
Banach  space  B.  Assume 

1.  E{|(p,  x*)|>  < for  all  x*  c 8* 

11.  E{(p,  x*)}  Is  continuous  In  x*. 

Then  there  exists  a unique  m c B such  that 

p 

E{(p,  x*)>  - (m  , x*) 

p 

Since  E{(p,x*)}  Is  a continuous  linear  functional  on  B*,  so 

of  B**  • B«  m fs  termed  the  mean  of  the  random  variable  p. 
p 

1.  m . « m + m , 

p+ir  p 

«.  Ihpll  < Ei||p||}. 

ill.  If  L:  3 ♦ B,  Is  bounded  and  linear,  then 

m,  ■ L m . 

Ip  p 

As  In  most  stochastic  processes,  the  mean  Is  not  our  prime  concern.  In  the  sequel 

we  thus  assume  that  all  random  variables  have  zero  mean.  For  the  definition  of  the 
variance  operator,  we  have  to  assume  the  following. 

1.  E{ | (p,  x*)  (t,  y*)|>  < -,  for  all  x*.  y*  t B*. 

( b) 

11.  E((p,  x*)  (»,  y*)>  Is  continuous  In  x*  and  y*. 

It  Is  easy  to  show  that  condition  (b)  implies  condition  (a).  Furthermore,  we  have  the 
following  lemma  to  facilitate  the  definition  of  the  variance  operator. 


It  Is  represented  by  an  element 
It  has  the  following  properties 
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Lam.  A continuous  bilinear  functional,  (x|y),  on  a Banach  spaca  3 Is  also 
boundad  (l.t.  there  axlsts  M c R such  that 

!(x|y)|/(| |x| M |y| |)  < N,  for  all  x.  y c B). 

Now  If  we  fix  y*,  than  E{(o,  x*)  (*-,  y#)>  Is  a boundad  linear  functional  on  B* 
(hanca  an  element  of  8**«8).  Thl s Implies  that  thare  exists  a unique  py#  c B such  that 

(p . x*)  (»,  y*)}  • (p^,  x*)  for  all  x#  « 8*. 

If  we  now  define  a mapping  <3  : 8*  *-  By 

ptf 


%.  r ■ V- 

It  can  easily  be  verified  that  Q Is  linear.  Morever,  (3  Is  bounded,  since 

01f  P« 


I IQeir  y* 1 1 • I Ip„*I  I ■ sup  y): 


llx*| 

< sup  "ItrLL-LLri.l  • M„y*„. 
x*  ||x*|| 


The  operator  Q Is  termed  the  covariance  operator  of  the  random  variables  p and  *. 
pir 

Covariance  operators  satisfy  the  following  conditons: 

1.  Q([_p)  (Kir)  " L ^pir  L and  * ar*  bounded  and 

linear  operators  on  3. 

11  * Lit  Q;  « Q s then 

O 00 

9„„  ' % * «„  * 9.0  * 9,. 

Qp  Is  called  the  variance  operator  of  p. 

111.  Q,„  ■ Q *.  In  particular,  q*  • q . 
pir  irp  p P 

1v.  q Is  positive;  l.e.  (q,,  y*.  y*)  • E((p,  y*)2>  > 0. 

M □ 


(2)  252L^2LJL^S2S£P-lSSSLJ£222S^i2^2!LZiCii£iS2i 

As  mentioned  In  the  previous  paragraph,  a reflexive  3anach  space  valued  random 

variable  has  a variance  operator  q which  Is  positive  and  self-adjoint.  Then  Thm.III.2 

p 


and  Thm.III.4  com  Into  the  picture  and  wo  have  tha  following:  Thare  exists  a Hilbert 

resolution  space  (H,E)  and  a left- factorization  (nrlnl  phase)  K:  (H,E)-*(B,F)  such  that 

Q » R 1C*.  Moreover,  tha  RKRS,  (H  , E),  which  corresponds  to  (Hn,nE)  In  Thra.III.4  Is 

also  assured  to  exist.  Since  wa  have  R(Q  ) c H c B,  one  of  tha  natural  questions  to 

0 0 “ 

ask  Is  whether  tha  random  variable  takas  values  only  In  H , and.  If  it  does,  what  can  wa  say 

p 

about  tha  original  random  variable.  The  answer  to  tha  first  part  of  the  question  Is  no,  and  a 
counter- 'example  has  bean  constructed.^  However,  If  we  happen  to  have  the  random  variable  takli.o 
values  In  Hp,  we  would  have  the  following  properties.  First  define  8Q  • R(k)S,  where 

1C:  H -►  8.  Then  K"1*:  B0  H and  we  have  (K”1")*:  H * B*.  Consider 
E{(p.  *)p(p»  y)fl>.  for  x,  y In  Hp. 

E<(p.  x)p(p,  y)p>  • Ef (1C“L  o,  K"L  x)h(K-L  p,  K’L  y)H> 

• E{(p,  (K’L)*  IC'L  x)B  (p,  (K"L)*  K_L  y)8  }. 


Note:  It"**  x and  (K'u)*  y are  elements  of  8*  , i.e.  linear  functionals 

on  8q.  By  the  Hahn-8anach  theorem^1,  there  exist  x*jand  y*  In  B*  such  that 

x*lao  * (K‘L)*  K'L  x. 
and  y*|R  « (IC'L)*  K’L  y. 


(p.  (K”L)*  K"L  x)a  - (p,  x*)f 
o 

(p.  Oc‘L)*  ic’L  y)«  ■ (p.  y*), 


Therefore, 


£{(p.  x)p(p,  y)p}  • E{(p,  x*)a(p,  y*)g} 

• (Qp  x*.  y*)g 

«■  (KK*  x*.  y*)8  - (K*  x*, 
■ (K"6x,  K“6  y)u  * (x,  y) 


as  such  the  random  variable 


43 

p has  tha  Identity  operator  on  H as  Its  variance  operator.  The  whole  Idea  can  be 
explained  by  the  following  example: 

Example  3.  Let  r • 5/12  and  h - (1,  (l/2)p,  (1/3)r.  ...)  c^,  where  ^ ■ 

{real  number  sequence  {y,}  J f y< 3 Is  finite).  Let  x be  a zero-mean  random  variable 

1-1 

taking  values  In  the  real  number  set  R and  let  Its  variance  be  1.  The  p - xh  Is  a random 
variable  taking  values  In  4 3.  The  dual  space  of  Is  43/2  snd  *3  ^ reflexive. 

(1)  For  all  z c ^/Z* 

E{(p,z)>  • E{(xh,  z)>  • E(x(h,z)> 

• (h,  z)  E{x}  - 0. 

Hence  p Is  zero-mean. 

(11)  For  all  z,  «t.4;  2/2 ’ 

E{(p,  z)(p,w)}  • E{(xh,z)(xh,w)}  ■ E(x2(h,z)(h,w)> 

• (h,z)(h,w)  £{x2)  ■ (h,z)(h,w) 

(Qaz.  ")• 

Thus  Q z ■ (h,z)h. 
p 

(111)  Since  the  range  of  Q Is  1 -dimensional , the  possible  factor  space  Is  the 

p 

real  set  R.  Let  K*  ■ 4^  -►  R be  the  functional  (h. * ) . 
then  K-R  -►  4-^  ■ 4^  Is  defined  by  the  following  equation, 

(Kd.z)  ■ (d,  X*  z),  for  all  dcR  and  zc^^* 

Since 

(d,  X*  z)  • d X*  z ■ d(h,z)  ■ (dh,z), 

we  have 

Xd  ■ dh,  for  all  dcR. 

For  all  z c ^2*  **  h4v* 

XX*  z • X(h,z)  • (h,z)  h ■ Qpz. 

Hence,  Q ■ XX*  . _ 

P 
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(1v) 


Furthermore,  we  have 


R*  def  RdC^F*)  - fR 
(.{0} 


* for  t > 0 
, for  t < 1 , 


RdCf*)  - K*Ft  C^3/23 

- (K^F*  z | z c 
{(h,Ft  z) | z c ^ 


, for  t > 1 
, for  t < 1. 


Hence,  £c  on  R Is  the  step  function 

U(t-1)  • f1  • t > 1 
to  , t < 1 

(v)  a.  When  E*d  ■ 0, 

case  1 (t  < 1)  ■ F*  ■ 0,  so  F^d  • 0; 
case  2 (t  > 1)  ■ E*  ■ IR  • 1,  so  d » 0. 
Thus  Fl  K d » 0. 


b.  When  Fl  K d ■ 0, 

case  1 (t  « 1)  ■ E*  ■ 0,  so  E*  d 


0; 


case  2 (t  > 1)  ■ F*K  d ■ Ffc  dh  • d F*  h • 0 

However,  F*  h - (1,  (7/2)r.  (T/3)r,  ...  (7/1)r,  0,  ...) 

* 0,  where  1 11  It. 

Hence,  d ■ 0.  Thus,  Et  d • 0. 

From  a and  b,  K Is  a left-mini  phase. 


(vl)  Hg  ■ K [R]  • {dh  ( d c R> 

For  d^  h,  d2  h c Hq,  the  Inner  product  Is  defined  as 
(d-j  h,  d2  h)g  • d1  d2 

(vl 1 ) Since  o takes  values  In  Hq  only,  o can  be  considered  as  a zero-mean  random 
variable  over  Hq  with  the  Identity  mapping  on  Hq  as  Its  variance.  This 
statement  can  be  verified  easily. 


fata,  slnca  the  stochastic  character  of  tha  given  random  variable,  0,  was  derived 
from  tha  scalar  random  variable,  x,  it  Is  appropriate  that  0 can  be  characterized 

completely  In  terms  of  Its  representation  In  the  one  dimensional  RKRS,  H . 

0 

j 

7.  Scattering  Operator 

Classically  In  network  analysis,  network  variables,  such  as  voltage  and  current, 
are  assuned  to  be  In  Hilbert  space.  Although  the  scattering  variables  are  a very 
useful  tool  in  network  analysis,  the  significance  of  the  normalizing  Impedance 
used  In  their  derivation  Is  not  clear.  Situations  have  occurred  where  we  have  to 
assume  that  the  network  variables  are  defined  In  8anach  space.  If  the  scattering 
variables  are  to  be  well  defined  here,  the  function  of  the  normalizing  Impedance 
should  be  the  transformation  of  network  variables  defined  In  Banach  space  Into 
elements  of  a Hilbert  space.  Theoretically,  It  Is  much  easier  to  work  with  Hilbert 
space.  Therefore,  the  significance  of  the  normalizing  Impedance  lies  In  the  fact  that 
It  transforms  a problem  defined  In  Banach  space  Into  a Hilbert  space  problem.  In 
this  section,  we  will  extend  the  Idea  of  scattering  variables  to  networks  with  their 
voltage  and  current  variables  defined  In  Banach  spaces  with  the  help  of  the  factor- 
ization theorems  developed  In  Chapter  III. 

Thinking  of  the  Impedance  Z,  or  the  transfer  function,  as  an  operator  from  a 
current  space  to  a voltage  space,  the  power  v j Is  a scalar  quanity.  Here  V denotes 
voltage  and  I current,  and  the  power  equality  Implies  that  the  voltage  plays  the 
role  of  a linear  functional  operating  on  the  current.  Thus,  Z may  naturally  be 
viewed  as  a mapping  from  a current  space  to  Its  dual,  a voltage  space.  Similarly, 
an  adnlttance  assumes  a dual  role  mapping  voltage  to  current.  In  that  case,  "I"  Is 
a linear  functional  over  the  voltage  space.  What  model  could  be  better  than  a re* 
flexlve  Banach  space  to  fit  the  structure?  Therefore,  In  this  work,  the  current 
and  voltage  spaces  are  chosen  to  be  dual  reflexive  Banach  spaces. 

Unlike  the  classical  case,  the  normalizing  Impedance  (or  admittance)  operator  is 


not  invertible.  In  general.  Besides  being  positive,  causal,  linear,  and  bounded. 
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we,  however,  also  assume  that  the  normalizing  Impedance  Is  1-1.  For  this  class 
of  operator,  we  have  the  following  special  factorization  theorems. 

Thm.  V.l.  Let  Z:  (B,  F)  ♦ (B*,  F*),  where  B Is  a reflexive  Banach  space, 
be  positive,  causal,  1-1,  linear  and  bounded.  Let  M ■ 1/2  (Z  Z*);  then  M is  also 
positive  and  furthermore,  self-adjoint.  Then 

1.  There  exists  a Hilbert  resolution  space  (H,  £)  and  a left- factorization 
of  M,  Kq:  (H,  E)  <*  (B*,  F*),  such  that  Ko  Is  left-mini  phase. 

11.  There  exists  a Hilbert  resolution  space  (H,  £)  and  a right- factorization 
of  M,  T0:  (8,  F)  ♦ (H,  £),  such  that  T0  Is  right-mini phase. 

Proof:  The  existence  of  the  left-  and  right-factorization  follows  from 
Thm.  III. 2 and  Thm.  III. 3. 

Note  here  that  we  use  the  same  Hilbert  space  for  left-  and  right- factor- 
ization. This  can  be  justified  from  the  proofs  of  Thm. III. 2 and  Thm.III.3  In 
Appendix. 

Thm.  V.2.  Let  M be  defined  as  above.  If  K:  (H,  E)  *»  (3*.  F*),  Is  a linear 
bounded  operator  such  that 

1.  1C  Is  1-1.  and  causal,  . 

11.  M ■ K K*, 

A*  A* 

then  there  exists  a linear  bounded  operator  U:  (H,  E)  ♦ (H,  E),  such  that 

a.  X • Ko  U,  (Kgls  as  defined  In  Thm.Y.l.) 

b.  U Is  causal  and  unitary. 


I 
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r 


U. 


Proof:  (1)  For  ill  y t H such  that  Ky  c Ko[H]t 
define  Uy  ■ Ko’L  K y. 

For  b c 3*  mo  have 
KK*b  ■ KoKo*b  c KoCH]. 

Honct  U Is  doflnod  ovor  IC*[B]  which  Is  donst  In  H. 

For  all  y c K*[B],  thoro  exists  x c B such  that  y ■ K*  x and  wo  havo 
l|Uy|  |2h  • 1 1 no"*-  K y ||2H  • ||to-L  It  It*  x ||2H 

• | |Ko"l  ko  Ko*  x ||2H  - | ( Ko*  x ||2h 

■ (Ko*  x,  Ko*  x)H 

• (x,  Ko  Ko*  x)g 

• (x.  K (*  x)g 

■ (K*  x,  1C*  x)H 

• llfx  ll*H-  I |y  ll2H 

Thoreforo,  U Is  Isometric  on  K*[B].  U can  thus  be  extended  over  H 

Isometrlcally. 

(11)  Similarly,  define 

7 y » K“L  Ko  y over  Ko*[B]. 

7 Is  Isometric  and  can  be  extended  over  H Isometrlcally. 

(Ill)  For  all  It  c H,  there  exists  sequence  (xl)  In  B such  that 

Ko*  xl  - h. 

This  Implies 

7 (Ko*  xl)  * 7 h 
and 

U 7 (Ko*  xl)  - U 7 h 


However, 

7 (to*  xt)  • K"L  Ko  Ko*  XJ 

• K"L  K 1C*  x1 

• K*  xi 

Hence 

U 7 (Ko*  Xj)  - Ko’L  K 1C*  xi 
• Ko*L  Ko  Ko*  x1 
■ Ko*  x^  ♦ h. 

Thenfore,  U 7 h ■ h and  U Is  an  onto  mapping.  Hanca,  U Is  unitary. 
(1v)  For  all  h c H,  than  exists  sequence  (x^ > such  that  K*  x t * h. 

We  also  have, 

Ko  U K*  Xi  ■ Ko  Ko_L  K K*  xf  • Ko  Ko"L  Ko  Ko*  x< 

• Ko  Ko*  Xf  • K K*  xi  - K h, 
but 

Ko  U 1C*  xi  * Ko  U h. 

Hence 

Ko  U h • K h for  all  h c H,  l.e.  Ko  U ■ K. 

(v)  Since  Ko  U » K,  we  have 

(F£)*  Ko  E£  U ■ (F£)*  Ko  U,  (Ko  is  causal) 

• (f£)*  K 

• (F£)*  K E£,  (K  Is  causal) 

■ (F£)*  Ko  U E£ 

• (F£)*  Ko  E£  U E£,  (Ko  Is  causal) 

Hence, 

(F£)*  Ko  E£  [E£  U - E£  U E£]  - 0. 

This  Implies  that  E£  U - S£  U E£  * 0,  since  Ko  Is  le*t-mini phase. 


Therefore,  U is  causal. 
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Thm.V.3.  Let  N be  defined  as  above.  If  T « (B,  F)  * (H,  £)  Is  a linear 
bounded  operator  such  that  1.  T Is  causal  and  has  dense  range. 

11.  N • T*  T 

Then  there  exists  W • (H,  £)  ♦ (H,  |)  such  that 

a.  T • W To 

b.  W Is  causal  and  unitary. 

Proof:  (1)  Oeflne  W*  y • (To*)*1"  T*  y,  for  y c T [8]. 

By  similar  argument  as  In  Thm.V.2.,  B*  can  be  proved  to  be 
Isometric  on  T[B].  Since  T[B]  Is  dense  In  H,  B*  can  be  ex- 
tended to  H Isometrlcally.  Also  like  In  Thm.V.2,  B*  can  be 
proved  to  be  unitary. 

(11)  Also  by  slmllary  argument.  It  can  be  proved  that 
To*  B*  • T*.  Hence  T • B To 

(111)  (Ft)*  To*  Et  B*  - (Ft)*  To*  B*.  (To*  Is  anti -causal) 

• (Ft)*  T* 

■ (Ft)*  T*  (T*  Is  anti-causal) 

- (Ft)*  to*  B*  Et 

■ (Ft)*  To*  Et  B*  Et 
Hence  (Fft)*  To*  (Et  B*  - Et  B*  Et)  -0. 

This  Implies  that  Efc  B*  - Et  B*  E?  • Q,  since  To*  Is  left- 
maxi  phase. 

Therefore,  B*  Is  anti-causal,  l.e. 

B Is  causal. 

It  Is  trivial  to  show  that  If  we  have  a causal  and  unitary  operator  U on 
(H,  E)  then  Ko  U Is  causal  and  (Kll)(ICU)*  ■ M.  Similarly,  If  we  have  a causal 
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«nd  unitary  operator  W on  (H,  |) , than  U To  Is  causal  and  (W  To)*(W  To)  » M. 

Tha  significance  of  thasa  facts  Is  that  thay  enable  us  to  choose  causal 
and  unitary  operators  U 4 W as  desired  In  order  to  make  tha  factorization  satisfy 
soon  additional  requirements.  Unfortunately,  the  proof  of  existence  and  the 
construction  of  these  U 4 H is  currently  beyond  our  reach,  even  though  It  Is 
trivial  to  do  so  In  the  classical  case.10  With  this  In  mind,  let  us  now  con- 
sider the  following  network: 


* f — - ■■  ■ 

k 

Vi' 

a a 

a.  n-port  network  series  loaded 


b.  optimal  matching  situation 


Figure  1. 

In  the  figure,  we  have 

1.  Ia,  rB  • a reflexive  Banach  space 

2.  Va,  Vg,  Vf  « B* 

3.  ZQ,  Z^  ■ (B,  F)  ♦ (B*.  F*),  linear  and  bounded. 

Although  the  circuit  diagram  Is  as  simple  as  shown  in  Figure  1,  there  are 

certain  requirements  for  the  diagram  to  be  well  defined.  They  are 

1.  Z„  + Z * is  1-1 
o o 

11.  Z0  ♦ Zj_  is  1-1 
ill.  Vg  « R(Zq  * ZQ*)  (1  R(Z0  + ZJ 


fcL- 


From  the  diagram,  we  have 


V ■ (Z_  + Z *1  X.  • (1  > Z. ) I 
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Ikflnt  Vr  - Va  - V1  and  Ir  - -(1^  - 1^,  than  we  have  Vp  • ZQ  Ip. 

Oaflna  Ip  » S1  Ij,  where  S1  is  call  ad  tha  currant-basis  scattaring  oparator,  than 

rr  * r1  * ’ *1  - * V*"  <Zo  * V>  't 


JB  " ♦ Z^)  *"  (ZQ  + ZQ*),  where  1^  is  tha  identity  napping  on  currant 


spaca  8.  Now  let  K and  T be  the  factorizations  of  1/2* (ZQ+  ZQ?)  as  defined  In 
Thm.IV.2  and  Thm.IV.3,  i.e. 


M • 1/2*(Z0  ♦ Z0*)  ■ a*  • T*  T. 

Oafina  a ■ K*  1^,  b » T Ip  and  b ■ S a,  where  S is  the  so-called  sacttartng  operator. 
Then  I ■ S1  1^  Implies  that 

I-R  b - S1  (1C*)"R  a,  so  b - T S1  (K*)^  a.  Hence 
S • T S1  (K*)‘R  • T (K*)‘R  -2  T (Zq  + ZL)*1  K 
• C - 2 T Ya  R,  where  C - T (*•)**  and  Ya  • (ZQ  + ZL)“U. 

In  order  to  have  a causal  scattering  operator  S,  wa  need  a causal  C.  However, 

C • T (K*rR  is  not  causal  in  general.  By  Thm.IV.2  and  Thm.IV.3, 

C ■ VI  To  (KoT*  U,  where  K , T denote  tha  left-  and  right- factorization  re- 
pacetlvaly.  Therefore,  tha  requirement  for  tha  selection  of  U 4 « is  to  make  C 
causal. 

Similarly,  consider  the  following  network,  an  n-port  network  parallel  loaded  by 
an  n-port  network. 


a.  parallel  loaded  n-port  b.  optimal  matching  situation 


Figure  2. 
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In  Figure  2 we  have 

1.  I4,  Ig,  I|  c B » reflexive  Banach  space. 

11.  V Va  B* 

111.  Yl,  Yq  : (8*.  F*)  - (B,  F),  linear' and  bounded. 

As  before,  certain  requirements  are  needed  for  the  circuit  diagram 
to  make  sence.  They  are 

1.  Y„  + Y,  Is  1-1 

0 L 

11.  Y ♦ Y * is  1-1 

O 0 

*g«R(V  n R(Y0+V* 

With  the  help  of  the  circuit  diagram,  the  following  equations  can  be  easily 
verl fled. 


*g  • (vo  * V>  V1  ■ n, ♦ V v0 
- t,  . r0  vi  • * »,  V 
ir  - 1,).  vr  V v,  - v,. 


I-  * V 

r or 


def 


S'  V«. 


where  SV  Is 


the  so-called  voltage-basis 


scattering  operator. 

s'  • (\  * V‘L  «v  - V ■ -'b*  ♦ z«  t’fo  * V>-  "h,r« 

2.  ■ <\  * V‘L‘ 

def  def 

a • Q*  V^t  b • P Yp,  where  P,  Q are  the  factorizations  of 

1/3  <Y0  ♦ YQ*)  • P*  P ■ q Q*  as  mentioned  In  Thm.IV.2  and  Thin. IV. 3. 
def 

b • S a,  where  S Is  the  scattering  operator. 

S - -?  (0*)“R  ♦ 2 P Za  g • -0  ♦ 2 P Z4  g.  where  0 - P (q*}“*. 
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Similarly,  In  order  to  have  a causal  S,  we  must  have  a causal  0.  However, 

0 • P (Q*)-r  » W'  P0  (Q0*)“R  U',  where  P0 » Q0  Is  the  right-  and  left- factorization 
of  1/2*(YQ  * YQ*).  Therefore  the  requirement  for  the  selection  of  W 4 U'  Is  to 
make  0 causal. 

One  of  the  most  useful  properties  of  scattering  variables  Is  that  they  give  a 
measure  of  the  optimal  transducer  power  gain.  To  see  that  this  property  still 
holds  for  our  generalized  scattering  operators,  let  us  consider  the  "power"  entering 
the  load  network.  For  the  series  loaded  network, 

ra  * -*r  * J1  * * T**  b + O'*)"*  a 

• Vr  ♦ *i  ■ zo  T‘Rb  * V (**>’***• 

So  the  power  entering  the  load  Is  given  by: 


(Of)'* 

a - T**1  b. 

Zo 

T^  b 

♦ V 

' (If)-*  «)8 

(Cic-)** 

a.  Z0*  (1C*) 

-R 

a)B  - 

(T^ 

b,  la  T*  b), 

♦((if)*1 

Ra.ZorR 

- (T* 

•"  6. 

V (**>'"  *)| 

(*.  k*l 

z0*  (1C*)-R 

a)H 

- (b, 

, (T*) 

'L  *0  T'R  b>H 

♦ (*. «' 

a Z0  T" R b) 

H " 

0c*L 

Zor 

R K 1 

b.  .)H 

- l/2(a,  IC“LZ0*  0C*fR  a)H  + 1/2(1C"L  Z0*  (K*)'R  a,  a)H 
- l/2(b,  (T*)’L  Z0  r R b)H  - 1/2((T*)'1  ZQ  T"R  b.  b)H 

• (a.  1/2*K“ L (Zq*  ♦ Z0)  (IC*)“R  a)H 

• - (b,  1/2‘CT*)^  (ZQ*  «•  ZQ)  T“*  b)H 

• (a,  K‘L  IC  K*  (C*)-R  a)H  - (b,  (T*)"1  T*  T T"*  b)H 

• (a.  a)H  - (b,  b)H  " 

• (a,  a)H  - (S  a.  S a)H 
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4.  Since  we  have 

Ip2  * ip  * ^ 

'•  V ! *i  * V * V 

For  all  x c £p,  y c £g,  we 

(y.  ip2  *)2  ■ (y»  *)? 


have 


Hence, 


Similarly, 


Therefore, 


I *1 

1-1  1 

(*.  y)« 


Ip2*  y ■ y.  for  all  y c l.e. 

r * ■ i 
xp2  12q* 


r * « r 

*2q  *p2* 


rpq  " *2q  *p2  " *2q  *2q*  " Ip2*  *p2* 


5.  It  can  easily  be  proven  that  Ip2  and  *2q  ar€  causa*  and  anti -causal  *n 
the  usual  time  structure  of  the  natural  numbers.  We  also  have 

M $8*  i (lnn  * I ) - I„ 

2 pq  pq'  pq 

» 

«T  * T aT  T* 

*p2  *p2  *2q  A2q  * 

6.  Using  the  formula  derived  for  scattering  operator  with  T • Ip2  and 


1C  ■ I2q»  we  have 


s ■ Ipj  <I2q*)'ft  * 2Ip2  (Ip,  * Il>‘L  Ij, 


• St  W*  * 2 St  (Ipq  * S’'"  2* 

* »*  * 2 JP2  «pq  * V*  ‘2q 

where  Is  the  Identity  mapping  on  L j. 

The  transformation  Involved  can  better  be  explained  by  the  following 
diagram: 


The  scattering  operator  Is  meaningful  when 

and 

*pq  * \ 1S 

Note  In  this  example  operator  I Is  not  bounded.  However,  boundedness 
has  nothing  to  do  with  the  derivation  of  the  scattering  operator.  Boundedness 
only  guarantees  the  existence  of  the  factorization.  Once  the  factorization 
Is  given,  the  derivation  for  the  scattering  operator  follows  through  accord- 
ingly. 

VI.  Conclusions 

The  research  herein  originated  with  a discussion  between  one  of  the  authors 
and  Professor  Harley  Flanders  concerning  the  underlying  mathematical  nature  of 
electric  networks.  It  was  observed  that  in  the  scattering  formalism  the 
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•rargy  dlsapatlon  of  a network  was  given  In  terms  of  the  norm  of  the  network 
variables  by 

e-  ii«ii2-  :ibii2 

whereas.  In  the  Irani  ttance  formalisms  the  energy  was  given  by  the  Inner 
product  (or  functional)  equality 


E - (v,l) 

Since  a network  Is  fundamentally  an  energy  processing  system,  one  might  Initially 
Interprete  these  equalities  as  Implying  that  the  scattering  variables  are 
naturally  defined  In  a Banach  space  since  only  the  norm  Is  required  to  define 
their  energy,  whereas,  the  Immlttance  variables  for  which  an  inner  product  Is 
required  to  define  energy,  are  naturally  defined  in  Hilbert  space.  In  fact, 
the  situation  Is  just  the  contrary.  The  Irani ttance  variables  may  naturally  be 
extended  to  Banach  space  by  working  simultaneously  with  B and  B*  whereas,  the 
development  of  the  present  paper  Indicates  that  the  scattering  variable  "live" 

In  a Hilbert  space  even  when  their  corresponding  Irani ttance  variables  are  de- 
fined in  8anach  space.  Indeed,  we  believe  that  the  primary  contribution  of 
the  present  work  is  the  observation  that  certain  problems  naturally  "live"  in 
Hilbert  space.  Moreover,  they  may  be  transformed  into  a Hilbert  space  even 
when  initially  defined  in  Banach  space. 


1.  Proof  of  Thm. III. 2: 

* 

A.  1.  By  Thm. III. 1,  there  exists  a Hilbert  space  H and  a 


A A 

linear  bounded  operator  K from  H to  3*  such  that 
» 

q • k k . 

II.  Define  H « R(ic’)  and  K « ic|H:  H * 8*.  then 
1C*:  B**-fl  ♦ H. 

111.  (K*  b,  x)H  ■ (b,  K x)3  » (b,  K x)&  ■ (K*  b,  x)H,  for  all 

b c B,  x t H. 

Since  K*  b,  K*  b t H,  so  K*  b • K*  b.  for  all  b c B. 
iv.  K K*  b ■ K K*  b ■ K K*  b « Q b,  for  all  b e B.  So 

k k*  « ic  k « q. 

v.  Since  H ■ K*[B]  ■ K*[B],  so  K*  has  dense  range.  So 
K Is  1-1. 

B.  1.  Oeflne  H*  ■ R(K*  Ft)  and  let  Et  be  the  orthogonal  projection 
on  H.  Then  (E*)2  « (E*)*  » E*. 

11.  Since  Hfc  becomes  R(1C*)  which  Is  H,  as  t * »,  so 
11m  E4  - IH  (Et  - IH  weakly) . 

III.  When  s <,  t, 

Hs  ■ R(IC*  Fs),  H*  » R(K*  Fc)  and  R(FW)  c R(^) . 

Since  Ft  b « Fs  b + FCs.t)  b«  for  all  b e 8,  and 
for  b'  c R(FS),  b*  • Fs  b',  so  • 

Ft  b*  - F*  (F5  b')  + F[s.t)  (F*  b‘)  • F5  b'  - b' . 


I.e.  b'  c R( Ffc) . So  R(!C*  F5)  c R(K*  F*)  and 

Hs  * R(IC*  Fs)  c R(K*  Fc)  - H*.  So  E*  Es  • Es  E:  « E5. 
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1v.  With  E defined  for  all  (— ,t),  t c R,  E can  be  extended  to 
s(R)  uniquely  to  be  a spectral  measure,  t.e.  a resolution  of 
Identity. 

v.  Since  £*00  ■ H*  ■ R(K*  F*),  by  definition,  so  K*  Is  a rlght- 
maxl phase  (Oef.II.5).  So  K Is  a left-mini phase  (Thm.II.8). 

C.  1.  Let  K:  (H,  E)  ♦ (8*,  F*)  be  another  left-mini  phase  factorization  of  Q. 
11.  Oeflne  W on  R(5t  ) by  W (X*  b)  - K*  b.  M Is  well-defined,  for  If 
K b ■ K a,  then 

K K*  a ■ Q a • Hlf  a - £ R*  b ■ Q b • K K*  b.  But  K Is  1-1 , so 
K*  a « K*  b. 

111.  For  any  b c B, 

2 2 

||W  K b||H  • | |K*  bf|H  - (K*  b,  K*  b)H 

- (b,  K 1C*  b)g  • (b,  K IE*  b)B 

* (K*  b,  T b)J*  - I |K*  b||#  • 


so  « Is  Isometric  on  R(K*).  Since  R(K*)  Is  dense  In  Hi,  W can  be 

* 

Ismetrlcally  extended  to  H. 

iv.  For  all  z c R(K*),  z ■ K*  x,  for  some  x c B,  hence 
KWz-KWK*x«KK*x-KK*x»Kz.  So 
K W ■ K over  R(ic*).  So  K W ■ K over  H via  continuity  as 
U extended  to  H. 

v.  By  the  same  argument,  there  exists  V:  H - H,  V Isometric,  such  that 
K V - K.  So  K W V - K,  and  hence 
H « L,  since  K Is  1-1.  Similarly,  V W » Ij. 
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vl.  With  K W ■ K,  we  have 

(F*)*  K E*  W • (F^)*  K W ■ (F*)#  lit  • (F*)*  ?!  & 
- (F*)#  K W E*  ■ (F1)*  K E*  W E*. 


(K,  K are  causal).  So 

(FV  1C  (Et  W - E*  W E*)  z - 0 for  all  ztH,  hence 
E*  (E*  W - E*  W E*  ) z ■ 0,  (K  Is  1 eft-mini  phase) , 
l.e.  E*  W - E*  W E*,  so  W Is  causal . 

vll.  Similarly,  V Is  causal.  But  W • V*  which  Is  anti-causal  (Thin. II. 6), 
so  U Is  memoryless. 

2.  Proof  of  Thnt.III.3.  

Oeflne  T » K*,  then  T*  ■ 1C,  and  define  ■ R(T  F^).  Then 

the  rest  of  the  proof  follows  as  In  Thm.III.2. 
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ABSTRACT 

A generalization  of  the  classical  Nyquist  stability  criterion  to  nonlinear  and 
time- varying  systems  is  obtained  via  an  appropriate  homotopy  argument  in  the  space 
of  causal  invertible  (possibly  nonlinear)  operators.  Although  the  resulting 
stability  test  is  only  a sufficient  condition  in  its  most  general  form  it  reduces 
to  the  classical  necessary  and  sufficient  Nyquist  criterion  for  linear  time- 
variant  systems  characterized  by  a transfer  function  or  transfer  function  matrix. 

Although  appearently  abstract  the  homotopic  nature  of  the  proof  proves  to  be  quite 
transpearent  and,  as  such,  many  of  the  classical  sufficient  conditions  for  non- 
linear or  time-varying  systems  can  be  derived  from  the  generalized  Nyquist  cri- 
terion by  simply  constructing  a homotopy  (continous  deformation)  from  the  given 
system  to  a ststem  which  is  known  to  satisfy  the  generalized  Nyquist  criterion. 

This  is  illustrated  via  a simple  derivation  of  the  Circle  criterion  as  a corollary 
to  the  generalized  Nyquist  criterion. 

INTRODUCTION 

When  one  discusses  alternatives  in  multivariable  control  the  classical  debate 
between  the  advocates  of  frequency  and  time  domain  techniques  usually  comes  to  the 
fore.  The  former  Is  highly  intuitive  but  restricted  to  linear  time-invariant  sys- 
tems whereas  the  latter  is  amenable  to  efficient  computational  procedures  and  is 
readily  extendable  to  nonlinear  and  time-varying  systems.  A third  alternative  is 
the  operator  theoretic  approach  wherein  the  system  is  modeled  by  an  operator  on 
Hilbert  space.  In  the  view  of  the  author  such  an  approach  to  the  control  problem 
achieves  the  best  of  both  the  time  and  frequency  domain  techniques.  Since  the 
operator  theoretic  model  is  defined  in  the  time  domain  the  resultant  control  tech- 
niques often  hold  for  nonlinear  and  time-varying  systems.  On  the  other  hand,  oper- 
ator theoretic  techniques  are  formally  quite  similar  to  the  operational  calculus 
associated  with  the  frequency  domain.  As  such,  the  intuitive  character  of  fre- 
quency domain  control  theory  often  carries  over  to  the  operator  theoretic  approach. 


♦This  research  supported  in  part  by  AFOSR  Grant  74-2631. 
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The  purpose  of  the  present  paper  is  to  illustrate  the  potential  of  the  oper- 
ator theoretic  approach  to  multivariable  control  via  the  derivation  of  a general- 
ized Nyquist  criterion  which  is  applicable  to  nonlinear  and  time-varying  systems 
modeled  by  finite  gain  operators  on  Hilbert  space.  Although  only  a sufficient 
l condition,  in  general,  the  technique  reduces  to  the  classical  necessary  and 

1 2 

sifficient  Nyquist  criterion  for  linear  time-invariant  multivariable  systems  * 
and  it  appears  to  be  "tight"  in  the  general  case. 

3 

9 Although  the  derivation  holds  in  an  abstract  Hilbert  Resolution  Space  for 

the  sake  of  brevity  the  present  discussion  will  be  restricted  to  the  case  of 
systems  defined  on  the  space  l.!)  composed  of  n- vectors  of  square  integrable  funct- 
t ions.  For  this  space  we  define  the  norm 

1.  Ilfil2  - / f(q)tf(q)  dq 

R 

and  a family  of  truncation  operators  P .L^  — *•  l!J  by 
f ~ " * C 

2.  (P‘f)  (q)  = /f<l)  * i ‘ 

l.  0 q > t 

An  operator  T : » l!)  is  said  to  have  finite  gain*  if  there  exist  constants  M 

and  N such  that 

3.  I I Tf  I I 1 M | | f | | + N 

for  all  f in  L^.  In  some  sence  the  constant  M plays  the  role  of  a norm  for  the  non- 
linear operator  T.  Of  course,  for  linear  operators  M may  be  taken  to  be  the 
norm  of  T with  N * 0 in  which  case  T has  finite  gain  if  and  only  if  it  is  bounded. 

In  the  nonlinear  case  if  an  operator  has  a finite  Lischitz  constant  then  it  is  al- 
so a finite  gain  operator  though  many  finite  gain  operators  do  not  admit  Lip- 

3 

schitz  constants. 

We  say  that  an  operator  T:^  * is  causal"*  if 

4.  P1!  • P^TP1 
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for  all  t.  It  is  easily  shown  that  the  causal  operators  are  closed  under  oper- 
ator addition  and  multiplication  and  limits  taken  in  the  topology  defined  by  the 
gain  constants  M and  N.  Unfortunately,  they  are  not  closed  under  the  operation  of 
operator  inversion.  A classical  example  of  this  is  the  unit  delay  whose  inverse 
is  a predictor.  Indeed,  the  question  of  determining  whether  or  not  T-1  is  causal 
from  the  properties  of  T is  completely  equivalent  to  the  question  of  determining 
whether  or  not  a feedback  system  is  stable  from  the  properties  of  its  open  loop 
gain1 *3>4,5,6.  SUCf,f  the  following  discussion  of  the  generalized  Nyquist  cri- 
terion will  be  formulated  in  terms  of  the  problem  of  determining  whether  or  not  the 
inverse  of  a causal  operator  is  causal,  the  solution  to  the  feedback  system  sta- 
bility problem  being  obtained  by  applying  these  results  to  the  return  difference 
operator.3 

THE  NYQUIST  CRITERION 

The  classical  Nyquist  criterion  is  usually  formulated  in  terms  of  the  degree 

of  the  system  frequency  response.  For  such  frequency  responses,  however,  their 

7 8 

degree  is  simply  a representation  of  their  homotopy  class  ’ and  hence  we  formu- 
late the  present  discussion  in  terms  of  homotopic  operators.  In  particular,  we 
say  that  operators  TQ  and  T^  are  homo topic  in  the  space  of  causal  invertible 
operators*,  C(0),  if  there  exists  a continous  operator  valued  function  T:I  C(0) 
mapping  the  interval  [0,1]  to  the  group  of  causal  invertible  operators  such 
that  T(0)  = Tq  and  T(l)  = T^ . Our  main  theorem  now  may  be  stated  as: 

THEOREM:  Let  Tq  and  T1  be  fininte  gain  operators  which  are 
homotopic  in  C(0).  Then  if  TQ  has  a causal  inverse  so  does  T^. 

The  proof  is  based  on  the  following  lemma  usually  known  as  the  small  gain  the- 
orem1,3. 


*Both  the  operators  and  their  inverses  are  assumed  to  be  finite  gain  but  the 
inverses  need  not  be  causal. 


Lemma:  Let  T be  finite  qain  causal  operator  for  which  M < 1.  Then,  if  the  oDer- 

I 

ator  (1  + T)  has  a finite  gain  inverse,  the  inverse  is  causal. 

I $ 

A proof  of  the  lemma  appears  in  reference  1 and  will  not  be  repeated  here. 


Proof  of  the  Theorem:  Let  T be  a homotopy  from  Tq  to  T^  and  assume  that  T ( tQ ) has 

a causal  inverse.  Now  let  It  - tl  < e and  write 

' o' 

5.  T(t)  * T(t0)  + (T(t)-T(t0))  - [1  + (T(t)-T(t0))T(t0)-1]  T( tQ) 

By  hypothesis  (T(t)-T( tQ) )T( tQ) is  causal  and  has  a gain  constant  M < 1 if  e is 
f chosen  sufficiently  small  (by  continuity  and  the  fact  that  T( t ) if  finite  gain). 
Moreover, 

6.  [1  + (T(t)-T(t0))T(t0)’1]  -1  * T(t0)T(t)_1 

exists  and  is  finite  gain  since  T(tQ)  Is  finite  gain  and  T(t)  has  a finite  gain  in- 
verse by  hypothesis.  As  such,  the  small  gain  theorem  implies  that 

7.  T(t)'1  « T(t0)-]  [1  + (T(t)-T(t0))T(t0)-1]  _1 

is  causal  since  it  is  the  product  of  two  causal  operators.  Finally,  since  the 
[0,1]  interval  is  a compact  set  one  can  piece  together  finitely  many  e -intervals 
to  show  that  T(l)~^  = Tj^  is  causal  if  T(O)"1  = T”^  is  causal,  thereby  completeing 
the  proof. 

APPLICATIONS 

Intuitively  the  theorem  states  that  the  property  of  a finite  gain  causal 
operator  having  a finite  gain  causal  Inverse  is  an  invariant  of  the  arcwise 
connected  component  of  C( 0)  in  which  the  operator  lies.  To  obtain  a test  for 
causal  invertibility  it  therefore  suffices  to  show  that  a given  operator  lies  in 
the  same  arcwise  connected  component  of  C( 0)  as  an  operator  which  is  known  to  ad- 
mit a causal  inverse.  In  particular,  the  following  corollary  reduces  to  the 
classical  Nyquist  condition  for  linear  time  inveriant  multivariable  systems. 
Corrollary:  Let  T be  finite  gain  operator  which  is  homotopic  to  the  identity 
operator  in  C ( 0) . Then  T has  a causal  inverse. 


Of  course,  the  identity  in  the  above  corollary  could  equally  well  be  replaced 

by  any  operator  satisfying  one  of  the  classical  sufficient  conditions  for  causal 

13  13 

invertibil ity  ’ ; operators  satisfying  the  conditions  of  the  small  gain  theorem 

3 3 7 

monotonic  operators  , operators  of  the  form  1 + S where  S is  strictly  causal  , 

etc.  Interestingly,  however,  each  of  these  classes  of  operators  lie  in  the 
same  arcwise  connected  component  as  the  identity  and  the  fact  that  they  admit 
causal  inverses  is  most  easily  derived  from  the  above  corollary  rather  than  con- 
versely. Another  such  class  of  operators  which  fall  in  the  same  arcwise  connected 
component  as  the  identity  are  the  causal  operators  for  which  zero  lies  in  the  un- 
bounded component  of  their  resolvant  set.  Indeed,  the  proof  that  such  operators 
admit  causal  inverses  appearing  in  reference  9 is  almost  identical  to  the  proof 
of  the  present  theorem  but  with  a restricted  class  of  homotopy.  In  fact,  the  pre- 
sent theorem  is  a simple  extension  of  the  earlier  result  though  considerably 
tighter.  In  particular,  the  result  of  reference  9 is  not  necessary  and  sufficient 
in  the  linear  time-invariant  case  and  assumes  that  the  operators  involved  admit 
finite  Lipschitz  constants. 

An  alternative  way  of  looking  at  the  above  theorem  is  as  a perturbation 
theorem  wherein  large  purturbations  are  allowed  so  long  as  they  are  conti nous  re- 
lative to  the  operator  topology  defined  by  the  operator  gain.  As  such,  it  should 
not  be  surprising  that  many  of  the  small  perturbation  results  of  classical  sta- 
bility theory  follow  from  the  generalized  Nyquist  criterion.  For  instance,  one 
may  derive  the  circle  criterion^  via  a two  step  homotopy.  First,  one  deforms  the 
nonlinear  term  to  a linear  (lying  in  the  middle  of  the  sector  associated  with  the 
nonlinearity)  and  then  one  deforms  the  resultant  linear  system  into  the  identity 
operator  via  the  classical  Nyquist  criterion.  Here,  a combination  of  the  sectoral 
bound^  and  the  requirement  that  the  spectrum  of  the  linear  part  of  the  system  lie 

outside  of  an  appropriate  disk  suffices  to  assure  that  the  operator  lies  in  C ( 0)  at 
every  point  in  the  homotopy  and  hence  justifies  the  application  of  the  Theorem. 


Of  course,  once  this  homoptopic  point  of  view  is  adopted,  numerous  generalizations 
become  appearent. 


CONCLUSIONS 

The  purpose  of  this  short  paper  has  been  two- fold.  First,  we  believe  that  the 
generalization  of  the  Nyguist  criterion  presented  may  prove  to  be  an  extremely 
powerful  tool  of  stability  theory.  Indeed,  we  conjecture  that  this  single  ele- 
mentary results  subsumes  most,  if  not  all,  of  classical  stability  theory.  Second- 
ly, however,  we  believe  that  it  illustrates  the  power  of  operator  theoretic  tech- 
niques in  control  which  have  the  potential  of  achieving  the  best  of  both  the  time 
and  frequency  domain  worlds.  Indeed,  such  techniques  yield  natural  and  intuitive 
generalizations  of  the  classical  frequency  domain  concepts  without  the  linearity 
and  time- invariance  restrictions  usually  associated  therewith. 
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ABSTRACT 

Classically  the  study  of  closed  loop  system  stability  is  approached  through  fre- 
quency domain  techniques,  e.g.  the  Nyquist  and  Hurwitz  criteria.  In  the  nonlinear 
case  frequency  response  is  not  well  defined;  however,  one  of  the  authors  has  re- 
cently shown  that  the  spectrum  of  a nonlinear  operator  can  be  used  in  lieu  of  the 
usual  Nyquist  plot  as  a means  of  generalizing  the  Nyquist  criteria  to  the  nonlinear 
case. 

Through  some  perturbation  techniques  we  characterize  in  this  paper  the  stability  of 
nonlinear  operators  by  the  more  accessible  "approximate  point  spectrum"  as 
opposed  to  the  entire  spectrum. 


I.  INTRODUCTION 

Recently  one  of  the  authors  demonstrated  that  the  stability  of  a closed  loop 
system  rests  squarely  on  knowledge  of  the  spectrum  of  the  operator  which  represents 
the  open  loop  gain.  (1)  The  system  may  be  nonlinear,  multivariable,  and/or  time- 
varying.  For  a linear  operator  representing  the  open  loop  gain  the  spectrum  and 
frequency  response  coincide,  however,  computation  of  the  spectrum  of  a nonlinear 
operator  is  not,  in  general,  a trivial  exercise,  (1)  (2) 

This  paper  shows  that  knowledge  of  the  more  easily  computed  approximate  point 
spectrum  is  adequate  to  answer  the  stability  question.  (3)  In  reference  1,  it  is 
shown  that  if  the  spectrum  of  the  operator  (representing  the  open  loop  gain  in  a 
unity  feedback  system)  does  not  encircle  "-1"  then  the  system  is  stable.  Essentially, 
this  is  equivalent  to  the  requirement  that  the  component  of  the  resol vant  (the 
complement  of  the  spectrum)  contain  the  point  "-1" — i.e.,  "-1"  is  not  disconnected 
from  infinity  by  the  spectrum.  It  is  shown  here  that  the  infinite  component  of  the 
resolvant  and  the  infinite  component  of  the  complement  of  the  approximate  point 
spectrum  are  identical.  Thus  if  the  approximate  point  spectrum  of  the  afore- 
mentioned operator  does  not  encircle  "-1",  then  the  system  is  stable. 

Finally  the  set  of  complex  numbers  {X  3 y(u))/x(ui) } , where  xU)  and  y(u>)  are  the 
Fourier  transforms  of  the  input  and  output  respectively,  is  shown  to  contain  the 
approximate  point  spectrum. 

*This  research  supported  in  part  by  Air  Force  Office  of  Scientific  Research  Grant 
AFOSR  74-2631. 
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In  some  cases  this  set  Is  very  large  and  can  be  the  whole  complex  plane.  This 
objection  is  offset  somewhat  by  the  fact  that  the  spectrum  of  trivial  nonlinear 
operators,  such  as  a squarer,  may  also  be  very  large. 


II.  THE  APPROXIMATE  POINT  SPECTRUM  AND  STABILITY 


All  operators  map  l!)  to  itself  unless  otherwise  specified.  An  operator  W is  causal 
If  whenever 

x(t)  * y ( t ) t < T;  f.g  in  (1) 

for  some  T,  then 


(Wx)(t)  - (Wg) (t)  t < T.  (2) 

The  norm  of  W,  ||W||,  is  the  usual  Lipschitz  norm.  W is  stable  if  it  is  both 
causal  and  bounded.  (4)  (5)  The  spectrum  of  an  operator,  W (possibly  nonlinear), 
is  the  set  of  complex  numbers  x,  such  that  the  operator  (X  - W)  does  not  have  a 
bounded  inverse.*  o(W)  denotes  the  spectrum  of  the  operator  W.  o(W)  is  a 
compact  set.  The  resol vant  set  of  W is  p(W)  which  is  the  complement  of  the 
spectrum  of  W in  the  complex  plane.  Clearly  p(W)  is  open. 


The  following  is  a recent  theorem  by  one  of  the  authors.  (1) 


THEOREM  1 ; Let  the  open  loop  gain  of  a (possibly)  nonlinear  feedback  system  be 
represented  by  a stable  unbiased  transformation  X,  mapping  L^  to  itself.  Then  if 
the  spectrum  of  K in  the  algebra  of  Lipschitz  continuous  unbiased  operators  does 
not  encircle  the  point  "-1",  the  feedback  system  is  stable. 


The  theorem  says  that  the  closed  loop  system  is  stable  if  "-1"  is  in  the  infinite 
component  of  the  resolvant.  For  the  case  of  a single  input  single  output,  linear, 
time  invariant  system  whose  open  loop  gain  is  characterized  by  the  frequency 
response  H(u>)  the  Nyquist  plot  for  H(ui)  is  precisely  the  spectrum  of  the  open  loop 
gain.  Hence  the  above  theorem  coincides  with  the  classical  Nyquist  test. 


Typically  the  spectrum  of  a nonlinear  operator  W is  difficult  to  compute.  A 
characterization  of  stability  using  the  approximate  point  spectrum  offers  a more 
accessible  route,  at  least  theoretically.  To  this  end  we  denote  the  approximate 
point  spectrum  of  the  (possibly)  nonlinear  operator  W as  n(W).  w ( W)  is  the  set 
of  all  complex  numbers,  x,  such  that  for  all  c > 0,  there  exists  an  x f 0,  such 
that  | |(x  - W)x|  ( £c  1 1 x ( | . it(W)  is  a closed  set  and  contains  the  point 
spectrum  (the  set  of  complex  X,  such  that  there  exists  x satisfying  Wx  = \x). 

We  denote  the  complement  of  ir(W)  by  y(W).  A complex  number,  y»  is  in  >(W)  if 


*The  symbol  (x  - W)  where  x is  a scalar  is  used  to  denote  the  operator  (xl  - W) 
where  I is  the  Identity  operator. 
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there  exists  an  e > 0,  such  that  1 1 (x  - W)x||  > e| |x| | , for  all  x f 0.  In 
proving  theorems  this  definition  seems  to  have  more  utility  as  in  the  following 
proposition. 

Proposition:  ir(W)e  cr(W). 

Proof:  Suppose  x t o(W).  We  show  X t ir(W)  implying  ir(W)  c o(W).  Since  x i o(W) 
we  have  (X  - W)"^  exists  and  is  bounded.  With  this  fact  consider  the  norm  of  x. 

INI  • ||(X  - W)-1(x  - W)  x 1 1 < ||(X  - W)-1 1 1 1 1 (X  - W)x|  | (3) 

Setting  e 3 1/ | | (X  - W)”^||,  we  conclude  that 

I ! (x  -W ) x | | > c||x||.  (4) 

Thus  X i ir(W)  as  was  to  be  shown. 

Since  ir(W)C  o(W)  we  have  p(W)c  Y(W).  y(W)  and  p(W)  are  open  sets  since  they 
are  the  complements  of  closed  sets.  As  with  any  set,  both  are  the  union  of  their 
connected  components.  Both  contain  a unique  infinite  component.  Necessarily,  the 
infinite  component  of  y(W)  contains  the  infinite  component  of  p( W).  A corollary 
to  the  following  lemma  shows  that  these  infinite  components  are,  in  fact,  identical. 

Lemma  1:  Let  W be  a (possibly)  nonlinear  operator.  Let  y.  be  a connected 

- — a 


component  of  y(W).  If  y contains  a point  in  p(W),  then  y c p (W ) . 

a a 

Proof:  Essentially  we  show  that  each  connected  component  of  p(W)  coincides  with 

a connected  component  of  y(W).  Let  y be  a connected  component  of  y(W).  Suppose 

a point  p is  an  element  of  both  ya  and  p(W).  Suppose  further  that  q is  any  other 

3 

point  in  y,  and  that  i is  a path  connecting  the  points  p and  q. 

a 

Since  p is  in  the  resol vant,  (p  - W)”^  exists  and  is  bounded.  The  task  is  to  show 
1 

that  (q  - W)~  exists  implying  y c p(W).  Combining  this  fact  with  the  above  pro- 
position,  we  will  have  every  connected  component  of  p(W)  coinciding  with  some 
connected  component  of  y(W). 

The  idea  of  the  proof  is  to  use  the  definition  of  y(W)  and  the  compactness  of  the 
path,  i , to  find  an  e-ball  about  the  point  p,  such  that  for  any  x in  the  e-ball, 

(X  - W)’1  exists  and  is  bounded.  It  turns  out  that  the  e-ball  depends  only  on  a 
single  constant.  Thus  a finite  number  of  e-balls  can  be  pieced  together  along 
the  path,  I,  so  that  the  arbitrary  point  q is  in  the  resolvant.  The  details  now 
follow. 

I 
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By  definition,  for  any  A in  y . there  exists  an  m,  >0,  such  that 

u A 

||(A  - W)  x||  > mx | | x | | (5) 

for  all  x in  l!J.  We  can  choose  mx  continuously  by  taking 

inx  a suptm^}  (6) 

where  the  sup  is  taken  over  all  mx  satisfying  the  above  inequality.  This  modifies 
the  inequality  to 

1 1 (X  - U)x||  > mx 1 1 x 1 1 . (7) 

Since  the  path  t Is  compact,  in,  achieves  its  minimum  for  some  A in  i.  Define  m * 

_ A 0 

mx  > 0.  Thus  for  all  A in  i we  conclude 

| | (A  - W)x| I > m| |x| I . (8) 

Consequently  if  A is  in  i and  (a  - W)"^  exists,  then 

| | (A  - W)"1 | | < 1/m.  (9) 

It  remains  to  show  that  (a  - W)”^  exists  for  all  A in  a.  Define  Sm(p)  = 

(A/|a  - p|  < m>  to  be  an  m-ball  about  the  point  p.  Let  A be  in  Sm(p),  then 

(A  - W)  » ((A  - p)  + (p  - W}}  » ((A  - p ) (p  - W)’1  + l)(p  - W).  (10) 

This  factorization  is  valid  since  (p  - W)“^  exists  and  ;s  bounded.  The  norm 

| | (A  - p)(p  - W)"1 | | < |A  - p|  | | (p  - W)_1 | | < m(l/m)  * 1.  (11) 

By  the  contraction  mapping  theorem 

C(A  - p)(p  - W)’1  ♦ l]'1 

exists  and  is  bounded.  (6)  Since  (p  - W)'1  exists  and  is  bounded,  the  same  is 
true  of  (A  - W)"1  for  all  A in  Sm(p). 

All  p in  i,  such  that  (p  - W)"^  exists  and  is  bounded  can  be  enclosed  by  an 

0 0 

e-ball,  S (p  ),  where  e depends  only  on  m,  which  depends  on  the  path  t.  Piecing 

c 0 -1 
a finite  number  of  e-balls  together  we  conclude  (q  - W)  exists  and  is  bounded. 

Since  y(W)  is  open,  each  component  of  y(W)  is  open.  Thus  for  any  point  q in 

there  exists  a path  and  a number  m * m(t),  such  that  every  point  in  y.  is  in  p(W). 

u 

Thus  the  theorem  is  proved. 

Corollary  1:  The  Infinite  components  of  p(W)  and  y(W)  are  identical. 
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Corollary  2:  p(W)  is  the  union  of  some  subset  of  the  set  of  connected  components 
of  y(W). 

Proof:  Pick  all  components  of  y(W)  which  contain  a point  in  p(W).  By  the  above 
lemma,  their  components  are  contained  in  p(W).  By  the  previous  proposition  the 
result  follows. 

Corollary  3:  Let  W be  a bounded  operator.  If  y(W)  has  only  one  component,  then 

a (W)  * 7T ( W ) • 

Proof:  Since  W is  bounded  y(W)  and  p(W)  contain  identical  infinite  components. 
Therefore  since  there  is  only  one  component  of  y(W),  p(W)  = y(W)  implying  that 
a(W)  ■ ir(W). 

With  these  statements  as  a base,  a modified  general  Nyquist  criteria  follows. 

THEOREM  2:  Let  the  open  loop  gain  of  a (possibly)  nonlinear  feedback  system  be 
represented  by  a stable  unbiased  transformation,  W,  mapping  L^  to  itself.  Then 
if  the  approximate  point  spectrum  of  K in  the  algebra  of  Lipschitz  continuous 
unbiased  operators  does  not  encircle  the  point  "-1",  the  feedback  system  is  stable. 

Proof:  We  remark  that  "not  encircle  -1"  is  equivalent  to  "-1"  in  the  infinite 
component  of  the  resol vant.  Thus  "-1"  is  in  the  infinite  component  of  y(W). 
Therefore  if  the  approximate  point  spectrum  does  not  encircle  "-1"  neither  does 
o(W)  and  conversely. 

III.  A COVERING  OF  THE  APPROXIMATE  POINT  SPECTRUM 
For  a linear  single  input  single  output  operator,  H,  the  frequency  response  H(w)  = 

A A A A 

y(o>)/x(tu) , where  y (oi)  and  x(ui)  are  the  Fourier  transforms  of  the  output  and  input 
functions  respectively,  is  the  spectrum.  In  the  nonlinear  case,  it  appears  that 

A A 

the  relevent  object  to  study  is  in  fact  {y(w)/x(u))>  since  it  offers  a covering  of 

A 

the  approximate  point  spectrum  even  though  H(a>)  is  undefined. 

Let  W be  a (possibly)  nonlinear  operator.  Suppose  Wx  = y.  Define  S(W)  = Closure 

- * A A A 

{ A | X 3 y(o»)/x(w)  where  y(u>)  and  x(u>)  are  Fourier  transforms}.  With  these  assump- 
tions we  have  the  following  lemma. 

Lemma  2:  ir(W)C  S(W). 

Proof:  Suppose  X l S(W),  then  since  S(W)  is  a closed  set 

|x  - ^-\  > e (12) 

X ( o» ) 

for  all  x f 0,  for  all  w,  and  some  e > 0.  Now  by  Parsevals  equality 


I 


♦ 


||(A  - U)x||2 
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IUx(«)  - y(u)  1 1 2 

(13) 

fUxM  - y(cu)  |2 

(14) 

fix  - *W|x*(u)|2d. 

J X(u>) 

(15) 

|X  - ^i|2  fU(u)|2d» 
x(o)0)  -«•' 

(16) 

2 

by  the  mean  value  theorem.  It  now  follows  that  (16)  is  greater  than  e ||x||.  The 
lemma  is  true. 

2 

Suppose  Wx  = x . Utilizing  the  global  inverse  function  theorem,  one  can  show  the 
spectrum  of  W to  be  the  whole  complex  plane.  (7)  (8)  Thus  one  could  expect  the 
set  S(W)  to  be  large.  For  appropriately  restricted  weakly  additive  operators  yU) 
is  well  defined  in  terms  of  x(u)  and  knowledge  of  the  operator  W. 

IV.  CONCLUSIONS 

Apparently  the  approximate  point  spectrum  is  the  interesting  object  of  study  in 
the  stability  question.  Moreover  knowledge  of  tt(W)  offers  a sufficient  condition 
for  when  a(W)  = ir(W).  Perhaps  a variation  will  offer  a necessary  and  sufficient 
condition.  Lastly  the  set  S(W)  covers  tt(W).  The  set  S(W)  is  intuitively  satisfying 
since  it  can  be  interpreted  as  the  frequency  gain  or  frequency  response  of  the 
system. 
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ABSTRACT 


The  usual  proof  of  the  Nyquist  Theorem  depends  heavily  on 
the  agrument  principle.  The  argument  supplies  unneeded 
information  in  that  it  counts  the  number  of  encirclements 
of  "-1".  Stability  of  a system  requires  an  encirclement 
or  a no-encirclement  test.  Using  homotopy  theory,  this 
paper  offers  a more  intuitive  approach.  We  believe  this 
approach  will  lead  to  practical  generalizations.  For  ex- 
ample, systems  characterized  by  several  complex  variables 
such  as  multi  deminsional  digital  filters. 


I.  INTRODUCTION 

This  paper  introduces  a characterization  of  the  Nyquist  criterion  using 
homotopy  theory,  a branch  of  algebraic  topology.  The  authors  emphasize 
the  intuition  and  motivation  for  this  approach.  The  hope  is  to  aid 
interested  readers  to  further  extend  and  apply  these  ideas.  In  this 
vein,  proofs  are  omitted  so  as  to  simplify  the  presentation.  Details 
can  be  found  in  the  references.  With  this  philosophy  in  mind,  let  us 
define  the  type  of  system  we  will  be  discussing. 

As  illustrated  in  Figure  1,  let  g(s)  be  a rational  function  in  the  com- 
plex variable  s (bounded  at  s = «**)  representing  the  open  loop  gain  of 
scalar  single  loop  feedback  system.  The  closed  loop  system  has  transfer 
function  h(s)  = g(s)/[l+g(s)].  The  closed  loop  system  is  stable  if  agd 
only  if  all  poles  of  h(s)  are  in  the  open  left  half  plane  denoted  by  Z_ 
(where  Z will  denote  the  entire  complex  plane). 

The  Nyquist  Criterion  states  that  the  closed  loop  system  is  stable  if  and 
only  if  the  Nyquist  plot  of  g(s)  (i.e.  the  image  of  the  Nyquist  contour 
under  the  map  g(*))  does  not  encircle  nor  pass  through  "-1".  If  the 
Nyquist  plot  passes  through  "-1"  there  is  a pole  on  the  imaginary  axis; 
if  the  Nyquist  plot  encircles  "-1" ,Qthere  is  a pole  in  the  open  right 
half  plane,  which  we  will  denote  by  Z+  (£+  will  denote  the  closed  right 

♦♦This  boundedness  condition  can  be  dispensed  with  4 is  added  only  to  ease 
the  exposition. 

♦Supported  in  part  by  AFOSR  Grant  74-2631 d 


81 


'1 


half  plane).  The  following  section  constructs  the  required  machinery 
of  homotopy  theory 


Figure  1 

II.  MATHEMATICAL  PRELIMINARIES  & BACKGROUND 
Basic  to  homotopy  theory  is  the  concept  of  a path.  A path  or  a curve  in 
the  complex  plane  is  a continuous  function  of  bounded  variation  (2)  y :[o,l]-7 
y is  a closed  path  if  y(0)  = yO)-  y is  a simple  closed  path  if  y is  a 
closed  path  and  has  no  self  intersections.  The  image  of  I * [0,1]  under  y is 
called  the  trace  of  y and  is  denoted  by  IyK 


Figure  2 

Two  closed  curves  yq  and  y-j  are  homotopic  in  Z if  there  exists  a continuous 
> function  rtlxl-^  such  that: 

Intuitively,  yq  is  homotopic  to  y-|  if  one  can  continuously  deform  y0  into 
y-| . Moreover,  it  is  easily  shown  that  the  homotopy  relation  is  an  equiva- 
lence relation.  (4)  (5) 

‘ Another  important  property  of  a closed  curve  is  its  index  or  degree.  The 

index  (2)  is  of  closed  curve,  with  respect  to  a point  "a"  not  in  { y > is: 

n(Yia)  3 27T  / (z~a)  ^ dz 

y 
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This  integral  measures  the  net  increase  in  angle  that  the  ray  r of  Figure  3 
accumulates  as  its  tip  traverses  the  trace  of  y. 


Figure  3 

Intuition  for  the  approach  stems  in  part  from  the  observation  that  n(y,-l)  = 0 
if  and  only  if  y is  homotopic  to  a point  in  Z-  { -1 } (cf.  prop.  5.4,  ref.  2). 

We  will  henceforth  refer  to  such  a y as  being  homotopically  trivial.  Conversely, 
y encircles  "-1"  if  and  only  if  y cannotbe  continuously  deformed  to  a point 
in  Z-  {-1}.  These  ideas  appear  to  indicate  that  the  Nyquist  encirclement 
condition  is  fundamentally  a homotopy  concept.  The  tuition  is  further 
reenforced  when  one  formulates  the  Nyquist  criterion  on  the  Riemann  surface 

A A 

(2)  (3)  associated  with  a map,  f(s).  Assuming  f(s)  is  analytic  on  Z¥  and 
bounded  at  s*®,  the  image  of  simply  connected  regions  in  Z+  are  simply 
connected  in  0.  To  illustrate  the  point,  let  Figure  4-a  be  the  image  of 
the  right  half  plane  under  f(s).  The  region  is  not  simply  connected. 

Figure  4-b  shows  the  "same  region"  as  it  might  appear  on  an  appropriate 
Riemann  surface.  Here  the  region  is  simply  connected. 


f 1 jure  4 

’•  i . -•  i ir ■»  r*t*»  NvTj’st  nlots  of 


» 
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A 

f(s)  In  2 and  on  the  Riemann  surface.  On  the  Riemann  surface  the  Nyquist 
test  becomes  an  obvious  triviality.  In  2 it  is  mathematically  more  delicate. 

Our  setting  uses  homotopy  theory,  a branch  of  algebraic  topology,  to  establish 
a topologically  invariant  relationship  between  a metric  space,  X,  and  an 
algebraic  group  called  the  fundamental  group  of  X,  denoted  by  w(X).  The 
relationship  is  topologically  invariant  in  that  homeomorphic  spaces  have 
isomorphic  fundamental  groups. 

Specifically,  the  fundamental  group  is  a set  of  equivalence  classes  of 
closed  curves.  Each  equivalence  class  consists  of  a set  of  curves  homotopical ly 
equivalent.  The  group  operation  is  "concatenation"  of  curves. 

For  example,  the  fundamental  group  of  2 consists  of  one  element,  i„,  the 
identity,  since  all  closed  curves  are  homotopic  to  zero.  If  X * 2 - {-11, 
then  't(X)  has  a countable  number  of  elements:  ix  (the  identity)  equal  to 
the  equivalence  class  of  all  closed  curves  not  encircling  "-1"  and  The 
remaining  elements,  un  (n  3 1 , 2,  3 . . . ) consisting  of  the  equivalence  class 
of  all  closed  curves  encircling  "-1",  n times.  Moreover,  concatenated 
with  uk  is  equal  to  the  element 

Now  let  X and  Y be  metric  spaces.  Let  f:  X-*-Y  be  locally  homeomorphic.  In 
particular,  assume  that  for  each  point  y in  Y there  exists  an  open  neighborhood 
G of  y such  that  each  connected  component  of  f"'(G)  is  homeomorphic  to 
G under  the  map  f.  Under  this  condition  X is  said  to  be  a covering  space 
of  Y.  (2)  (4)  Also  let  *(X)  and  t(Y)  be  the  fundamental  groups  associated 
with  X and  Y respectively.  With  these  assumptions,  f effects  a group 
isomorphism  (i.e.  a one  to  one  into  mapping  preserving  group  operations) 
between  it ( X ) and  a subgroup  of  n(Y)  as  in  the  following  diagram  (4)  (5) 


X 


1 1 < K) 


V 

/ 


-X 


M c.  \\  V?) 


P 

Figure  5 


F is  the  functor  which  establishes  the  relationship  between  a topological 

space  and  its  fundamental  group.  Finally  let  us  distinguish  between  a 

critical  point  and  a critical  value.  A point  zQ  in  l is  a critical 

point  of  a differentiable  function  f if  f'(zQ)  * 0.  A critical  value  of 

f is  any  point  w = f(zQ)  whenever  zQ  is  a critical  point. 

flow  suppose  f:!Z  + 2 is  a rational  function  whose  set  of  poles  is  P = 

{P^>  . ..,  Pnl.  Let  Q = (q-j , ....  q^}  be  the  set  of  all  points  in  % such 

that  f(q.)  is  a critical  value  of  f.  Note  that  there  may  be  q.'s  which 
' 2 1 

are  not  critical  points.  To  see  this  consider  g(z)  *z  (z-a).  g (0)  3 0 

implies  "0"  is  a critical  vlue  of  g,  but  g(a)  = 0 with  g‘(a)  f 0. 

Finally,  define  T * { t ^ 1 1 ^ - f”^(-l),  i-1,  ....  n}.  Note  also  that  since 

f is  a rational  function,  P,  Q and  T are  finite  sets.  Define  X » Z - 

{PuQuT}  and  define  Y = f(X). 

Lemma  1 : Under  the  above  hypothesis,  X is  a covering  space  of  Y.  This 
leads  to  the  following  corollary. 

Corollary:  The  fundamental  group  n(X)  of  X is  isomorphic  to  a subgroup 
N of  H(Y). 

This  corollary  says  that  a closed  curve  in  X is  homotopical ly  trivial. 

III.  THE  SCALAR  CASE 

Let  g(s)  be  as  described  in  the  introduction.  Appropriately  define  the 
sets  P,  Q,  and  T and  the  spaces  X and  Y so  that  X is  a covering  space 
of  Y.  Also  as  per  reference  (10)  and  Figure  6,  construct  the  ugly 
Nyquist  contour,  and  the  usual  Nyquist  contour,  r,  where  r:I  - Z u W 


Figure  6- (a ) Figure  6-(b) 
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X indicates  a point  of  P;  $ indicates  a point  of  Q 

Lemma  2:  Under  the  above  assumptions  on  g and  xR,  h's)  is  stable  if  and 
only  if  the  path  go\R  does  not  encircle  "-l”.  (10) 

At  this  point  we  must  establish  this  lemma's  connection  with  the  classical 
Nyquist  criterion.  To  this  end  we  compare  the  information  of  the  flyquist 
plot,  gor  with  the  "ugly"  Nyquist  plot,  goxR. 

Lemma  3:  Let  n be  the  number  of  poles  of  g in  £ , then 


irT ) (z"1)’ldz  3 hr  \ (z-ir1dz  +n 

gor  go\R 

These  three  lemmas  give  rise  to  the  following  theorem. 

Theorem  1 : Let  g(s)  be  as  above.  Then  h(s)  is  stable  if  and  only  if  the 
Nyquist  plot  of  g(s)  does  not  pass  through  "-1"  an^j  encircles  "-1"  exactly 
n times  where  n is  the  number  of  poles  of  g(s)  in  Z+. 

IV.  MATRIX  CASE 

Let  the  entries  of  an  nxn  matrix  G(s)  be  rational  functions  in  the  complex 
variable  s.  Suppose  G(s)  characterizes  the  open  loop  gain  of  the  single 
loop  feedback  system  of  Figure  7. 


Figure  7 

x(s)  and  y(s)  are  n vectors  whose  entries  are  also  rational  functions  of 
s which  represent  the  input  and  output  of  the  system  respectively. 

A * 

This  article  assumes  each  entry  of  G(s)  is  bounded  at  X = « . Thus  G(s) 
as  a mapping,  G(*):£  -*■  £nxn,  is  analytic  on  £ except  at  a finite  number  or 
points,  the  poles  of  its  entries. 

For  Figure  7 to  be  well  defined  we  require  that  det  [I+G(s)]*0  Thus  there 
exists  a closed  loop  convolution  operator,  H,  such  that  y = H*x.  More- 
over the  Laplace  transform  of  H,  H(s)  satisfied 

H(s)  = G(s)[I+G(s)]-1 


For  this  system  to  be  stable,  H(s)  must  have  all  its  poles  in  Z and 
have  all  its  entries  bounded  at  s = «. 

Under  the  assumptions  on  G(s),  the  following  factorization  is  valid: 

G(s)  = N(s)D-1 (s) 

where  N(s)  and  D(s)  are  right  co-prime,  polynomial  matrices  in  s with 
det[0(s)]^0.  Moreover  Sq  is  a pole  of  G(s)  if  and  only  if  it  is  a zero 
of  det[D(s)].  (9) 

Dosoer  and  Schulman  (3)  have  shown  that  the  close  loop  operator  H is 
stable  if  and  only  if  det[N(s)+D(s)]/*0  for  s in  and  det[I+G(-)]j*0. 

Using  this  fact,  we  state  and  prove  the  following: 

Theorem  2:  H is  stable  if  and  only  if  (1)  the  Nyquist  plot  of  det[N(s)+ 

D(s)]  does  not  encircle  nor  pass  through  "0",  and  (2)  det[I+G(  )]^0.  (10) 
Observe  that  if  one  assumes  the  open  loop  gain  to  be  stable  (i.e.  G(s) 
has  all  poles  in  Z+)  then  det[I+G(s)]  in  the  above  theorem.  This  follows 
since  for  all  s in  £+,  det[N(s)+0(s)]=  det[I+G(s)]  det[D(s)]  with  det[D(s)lfO. 
Thus  in  Z+  det[N(s)+D(s)]  has  a zero  .if  and  only  if  det[I+G(s)]  has  a zero. 
Finally,  it  is  worthwhile  to  point  out  the  relationship  between  the  above 
formulated  multivariable  Nyquist  criterion  and  that  formulated  by' Barman 
and  Katznelson.  For  this  purpose  we  let  X.(iw);  j=l n;  denote  the  n 

A J 

eigenvalues  of  g(iw).  In  general  parameterization  of  these  function  by  iw 
is  not  uniquely  determined  but  one  can  always  formulate  such  a function. 
Moreover  these  functions  are  piecewise  analytic  and  can  be  concatenated 
together  in  such  a way  as  to  form  a closed  curve  which  Barman  and  Katzel- 
son  term  the  Nyquist  plot  of  G(s). 

Now,  since 

n 

det[I  + G(iw)]  = It  [1  + x.  (iw)] 
j=l  J 

and  the  degrees  of  a product  is  the  sum  of  the  degrees  of  the  individual 
factors  and  also  equals  the  degree  of  the  cancatonation  of  the  factors,  the 
degree  of  the  Barman  and  Katznelson  plot  with  respect  to  "-1"  coincides 
with  the  degree  of  our  plot  with  respect  to  "0".  As  such,  even  though 
the  two  plots  are  different  their  degrees  coincide  and  hence  either  can 
be  used  for  a stability  test. 

Acknowledgement:  The  authors  would  like  to  acknowledge  the  contribution 
of  Or.  John  Murray  (Dept.,  of  Mathematics,  Texas  Tech  University)  whose 
continuous  flow  of  counter  examples  shaped  the  ideas  presented  herein. 
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R.  A.  DeCarlo,  R.  Saeks,  and  M.  Strauss 
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Lubbock,  Texas 


ABSTRACT 

Using  two  classic  theorems  (one  of  Mackey  and  another  of  Strone)  and  a 
recent  result  of  Masani  and  Rosenberg,  this  paper  pieces  together  a generalized 
frequency  response  theory  for  an  abstract  Uniform  Resolution  Space.  The 
present  theory  assimilates  past' work  as  done  by  Falb,  Freedman,  Anton,  Masani 
and  Rosenberg,  and  one  of  the  authors.  The  results  of  the  paper  are  not  new, 
but  are  merely  a rearrangement  of  subtleties  uncovered  by  the  aforementioned 
authors.  An  interesting  consequence  of  this  work  was  that  an  abstract  Uniform 
Resolution  Space  has  both  a "time  transform"  and  a "frequency  transform".  Such 
a duality  is  not  readily  identifiable  in  an  L^  function  space  since  the  time 
transform,  there,  is  the  identity. 

INTRODUCTION 

Fourier  analysis  is  basic  to  the  design  and  understanding  of  physical  systems. 
The  property  that  convolution  in  the  time  domain  maps  into  a product  in  the  fre- 
quency domain,  yields  a theory  both  practical  and  aesthetically  pleasing.  This  rote 
provides  what  is  hoped  to  be  a generalized  frequency  response  theory  for  arbitrary, 
closed,  linear,  time  invariant  operators  on  a uniform  resolution  space.  Previous 
attempts  at  providing  a general  frequency  theory  have  illuminated  numerous  sub- 
tleties, yet  still  appear  inadequate  for  one  reason  or  another.  Interestingly 
enough,  the  mathematics  necessary  for  such  a synthesis  is  well  entrenched  in  the 

literature.  This  paper  merely  pieces  these  results  together  and  reinterprets  them 

3 9 12  21 

in  light  of  the  work  of  Falb,  Freedman,  Masani,  Rosenberg  and  Saeks.  ’ ’ ’ 

This  research  supported  in  part  by  Air  Force  Office  of  Scientific  Research  Grant 
AFOSR  74-2631. 


f 


91 


Classical  Fourier  analysis  consists  essentially  of  two  fundamental  ideas— 
the  idea  of  a "transform"  from  time  to  frequency  and  the  property  of  a time-in- 
variant mapping  to  a product  of  functions  in  frequency.  We  desire  a Fourier 
representation  for  time  invariant  operators  defined  on  an  appropriate  space. 

Two  avenues  arise.  A traditional  approach  uses  a Fourier-like  integral  to  obtain 
the  representation..  In  an  abstract  approach,  the  Fourier  representation  is  a 
spectral  representation  of  the  abstract  operator  relative  to  an  appropriate  spec- 
tral measure.  This  road  is  both  more  general  and  eliminates  the  need  for  a 

specific  representation  of  the  operator. 

3 5 

Falb,  Freedman  and  Anton  * developed  a generalization  closely  paralleling 
the  classical  theory.  The  formulation  considers  Hilbert  space-valued  L- 
functions  (square  integrable  relative  to  the  Haar  measure),  defined  over  a 
locally  compact  abelian  (LCA)  group,  6,  and  operators  which  are  characterized 
by  an  Lj  convolutional  weighting  function.  The  theory  is  highly  representation- 
dependent  and  fits  awkwardly  into  the  setting  of  an  abstract  resolution  space. 

In  fact,  the  identity  and  unit  delay  are  not  admissable  to  the  theory.  The  major 
advantage  is  that  one  obtains  an  operator-valued  Fourier  representation. 

Masani  and  Rosenberg1^ ^ use  a spectral  theoretic  vehicle  to  alleviate  the 
difficulty  of  a specific  representation  of  the  operator.  Moreover,  the  theory 
settles  nicely  into  an  abstract  setting.  Yet,  the  frequency  response  is  always 
scalar-valued,  even  in  the  multivariable  case,  and  the  concept  of  a "transform" 
is  absent. 

21 

Finally,  Saeks  has  a Masani-like  development  whose  Fourier  representation 
assumes  values  in  a suitably  restricted  class  of  operators.  The  advantages  are  the 
compatibility  with  abstract  spaces  and  an  operator-valued  frequency  response.  Yet 
still,  the  concept  of  a transform  is  missing  and  major  existence  questions  are 
still  present. 
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The  structure  of  the  present  theory  rests  on  the  classic  theorems  of 
Mackey^  and  Stone*  and  a recent  theorem  of  Rosenberg  and  Masani^.  With  this 
comment,  we  define  the  setting. 

UNIFORM  RESOLUTION  SPACE 

A resolution  space  is  a pair,  (H,E)  where  H is  a Hilbert  space  and  E is  a 
spectral  measure  cn  an  ordered  LCA  topological  group,  G.  On  an  ordered  LCA  group, 
a spectral  measure  determines  a resolution  of  the  identity,  and  conversely.  Thus, 
it  is  advantageous  to  work  with  the  resol tuion  of  the  identity  Et  = E([-®,t]), 
rather  than  with  the  spectral  measure  E,  as  illustrated  at  the  end  of  this  section. 

As  an  example,  consider  the  Hilbert  space,  l_2,  together  with  the  truncation 
operator,  Et,  defined  as 

(£tx)(q)  * -fx(q)  q < t 
l 0 q > t 

or  equivalently,  the  spectral  measure,  defined  via 


(E(S)x)(q) 


x(q)  q e 3 
. 0 q i B 


for  all  Borel  sets  B. 

In  addition  l 2 admits  a group  U of  shift  operators  U*,  defined  as 


(l^xMq)  = x(q  - t). 

Thus,  the  concept  of  time  invariance  is  well  defined  in  a classical  1_2  setting. 

In  general,  a resolution  space  lacks  the  concept  of  time  invariance.  Such  a 
property  requires  an  extension  of  the  concept  of  the  L2  "time-shift".  A group  of 
such  operators,  in  general,  fails  to  exist  in  an  arbitrary  resolution  space. 

In  particular,  we  seek  a strongly  continuous  group  of  unitary  operators 
(i.e.,  Ut_s  * Ut(Us)"^  for  all  t and  s in  G),  such  that 


UtE(B)  = E(B  + tJU* 


« 
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for  all  t in  G and  Borel  sets  B.  A resoltuion  space,  together  with  such  a 
group  U of  shift  operators,  Ul,  is  a Uniform  Resolution  Space  (URS),  denoted 
by  the  triple  (H,E,U). 

Underlying  each  URS  is  an  ordered  LCA  group,  G,  which,  for  our  purposes,  is 
time.  Associated  with  G is  a "character  group",  G,  which  is  the  group  of  con- 
tinuous homomorphi sms  from  G into  the  multiplicative  group  of  complex  numbers  of 

A 

magnitude  one.  Note  that  G is,  in  general,  not  ordered. 

In  like  manner,  attached  to  each  URS  (e.g.,  (H,E,U)),  defined  over  G,  is  a 

A A AAA 

"dual"  character  space  (H,U,E)*,  defined  over  G.  E and  U are  a spectral  measure 
and  a group  of  shift  operators,  respectively,  defined  via  the  two  equalities 

ut  - fgfY.-tlffilr)  t c G 
and 

^ ■ Jg(Y.-t)dS(t)  Y c G. 

Here,  (y,-t)  denotes  the  complex  number  of  magnitude  one,  resulting  from  the  oper- 

A 

ation  of  the  character  y in  G acting  on  -t  in  G,  and  where  the  integral  is  the 
Lebesque  integral.  Stone's  theorem(4)  assures  the  existence  and  uniqueness  of  E 

A 

and  U. 

A A A 

Oddly,  the  character  space  (H,U,E)  is  not  a resoltuion  space  since  G is  not 

A A 

ordered.  However,  (H,U,E)  displays  all  the  resolution  space  properties  which  do 

A 

not  depend  on  the  ordering  of  G.  In  fact,  by  Stone's  theorem(6),(7),(12),  U is  a 
group  of  shift  operators  for  E,  satisfying  the  imprimitivity  equality  over  G--i.e., 

u^ls)  - £(3  * Y)ur. 

*We  have  adopted  the  ordering  (H,U,E)  because  via  Stone's  Jheorem,  E and  U contain 
the  same  information.  Moreover,  U and  E do.  Thus,  (H,U,E)  rather  than  (H,E,U). 
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For  our  purposes,  the  character  group  plays  the  role  of  frequency. 

Now,  the  physical  properties  of  causality,  memorylessness,  time  invariance, 
etcetra,  have  precise  descriptions  in  the  uniform  resolution  space  structure. 

In  particular,  for  bounded  operators,  T,  on  (H,E,U),  causality  is  equivalent  to 
E%  * EtTEt;  anticausality,  to  EtT  * E^TE^.*;  memorylessness»  t0  E%  = 

which,  in  turn,  is  equivalent  to  T,  being  both  causal  and  anticausal.  Since 
memorylessness  is  a symnetric  concept,  it  has  an  analog  in  the  character  space, 

A A 

(H,U,E),  whereas  causality  does  not.  Because  of  this,  we  say  a bounded  operator, 

T,  is  time  invariant  if  E(B)T  * TE(B)  for  all  Borel  sets  B of  G.  Via  Stone's 
theorem,  this  is  equivalent  to  11%  = TU*  for  all  t in  G.  Clearly,  we  emphasize 
the  character  space  in  the  definition  of  time  invariance. 

In  the  case  of  unbounded  operators,  (e.g.,  the  derivative  operator),  T is 
causal  if  E%  c E^TE***;  T is  anticausal  if  EtT  _ EtTEt;  T is  memoryless  if 
E%  c TE%  and,  finally,  although  somewhat  non-intuitively,  T is  time  invariant 
if  E(B)T  c TE(B)  for  all  Borel  sets  B in  G,  where,  again,  we  emphasize  the 
definition  in  the  character  space.  For  unbounded  operators.  Stone's  theorem,  in 
general,  does  not  yeild  an  equivalent  statement  (such  as  U%  = TU%  in  the  original 
resol tui on  space.  However,  for  the  case  of  linear,  single-valued,  closed  oper- 
ators with  domain  dense  in  H,  then  U%  « TU*  if  and  only  if  E%  ^ TEt.®’^  The 
fundamental  role  of  the  character  space  becomes  more  clear  in  the  following  section. 
EQUIVALENT  SPACES 

In  this  section,  Mackey's  theorem  verifies  an  equivalence  between  an  abstract 
URS,  (H.E.U)  and  a function  space,  (L2(G,K),  Xg.o1).  Now,  the  relevant  infor- 
mation contained  in  (H,E,U)  is  also  contained  in  (H,U,E).  Thus,  applying  Mackey's 
theorem  to  (H.E.U)  (under  the  guise  of  (H,U,E)),  another  equivalence  to 
(L2(G,K),  oy , Xg)  exists.  Furthermore,  (L2(G,K),  Xg,  at)  and  (L2(G,K),  <jy,  Xg)  have 
and  affinity  via  Stone's  theorem. 

*Et  - E((t,  -))  - I - E . 

**For  an  unbound*!  operator,  T,  onfca  resolution  space,  (H.E.U),  the  detain  of  E*r 

Is  smaller  than  the  domain  of  TE  . As  such,  the  containments  Indicate  that 

where  the  domains  of  Em  and  TE*  coincide,  then  E^T  - T£t.  * f 


A 
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After  numerous  references  throughout  this  development  to  the  above  authors, 
we,  at  last,  precisely  state  their  results.  Hopefully,  this  will  facilitate 
understanding  of  the  maps  between  the  various  spaces,  hinted  to  in  the  above 
paragraph.  The  following  is  a statement  of  Stone's  theorem  for  LCA  groups.^’ ^ 

A 

Suppose  G is  an  LCA  group  and  G,  its  "dual"  character  group;  let  (y,  -t)  be 
the  complex  number  of  magnitude  one,  resulting  from  the  operation  of  y in  G on  -t 
in  G;  define  e(e)  as  the  a-algebra  of  Borel  sets  of  G(G);  finally,  let  {Ut|t  in  G} 
({UY|y  in  G})  be  a strongly  continuous  group  of  unitary  operators  on  a complex 

A 

Hilbert  space,  H,  onto  H.  Then,  there  exists  a unique  spectral  measure,  E(0 
(E( - ) ) » for  H on  e(e) , such  that  all  t in  G (y  in  G), 


or 


§(y.  -t)dE(r) 


- f g(n  -t)ds(t). 


The  Initial  task,  now  is  to  construct  an  equivalence  between  two  L2  spaces  via 
this  theorem. 

Consider  the  URS,  (L2(G,K),  Xg,  o*),  and  the  character  space,  (L2(G,K), 
oY,  Xg),  where  G is  an  LCA  group;  G,  its  character  group,  Xg,  the  characteristic 
function  of  the  Borel  set  B in  E;  ot,  the  classical  shift  operator  (i.e., 
(atf)(q)  ■ f(q  - t)*,  and,  lastly,  the  measure  on  the  space  will  be  the  Haar 


measure,  m. 

A 

Now,  the  Fourier  transform  maps  L2(G,K)  to  L2 ( G,K)  in  such  a manner  that 
Xg  Is  taken  to  the  spectral  measure  on  L2(G,K),  whose  integral  is  oY.  Similarly, 

4 A 

a maps  to  the  unitary  group  on  L2(G,K),  whose  associated  spectral  measure  is  Xg. 
As  such,  (L2(G,K),  cy,  Xg)  is  the  Fourier  transform  of  the  character  space  for 
(L2(G,K),  Xg,  ol). 


{°t,e  ,n  c>  *.ndw^a,B  in  ,et  of  3cp•,  of  G)  serv»  as  the  strT>n51y 
continuous  group  of  shift  oparators  and  th*  spactral  mtasura.  r»sp#ct1v*1y. 
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Interpreting  this,  we  nave  Xg  completely  determining  c^and  Xg,  completely 
specifying  ot.  The  link  between  these  two  spaces  and  an  abstract  resolution 
space  is  Mackey's  theorem.  The  statement  of  his  theorem  for  LCA  groups  follows; 

Let  E be  a spectral  measure  on  the  Borel  sets  of  an  LCA  group,  J,  and  let 
U be  a strongly  continuous,  unitary  representation  of  J,  such  that 

E(B  + tJU*  - ll^EfB)** 

for  all  Borel  sets  B of  J and  all  t in  J;  then,  there  exists  a unique  Hilbert 
space,  K,  and  a unitary  transformation,  M, 

M:H— l2(j,k,b) 

such  that 

(1)  MEtBjM'1  - Xg  for  all  Borel  sets 
of  J;  and 

(2)  MliV1  - a 1 for  all  t in  J, 

where  K is  a complex  Hilbert  space  and  m,  the  Haar  measure. 

For  an  arbitrary  URS,  (H,E,U),  defined  over  an  LCA  group,  G,  by  design, 

E and  U satisfy  the  imprimitivity  equality.  Moreover,  the  character  space, 

- « - 20 
(H,U,E),  possess  the  property  that  E and  U satisfy  the  imprimitivity  equality. 

Finally,  it  is  clear  that  (L2(G,K),  Xg,  ot)  and  (L2(G,K),  cY,  Xg)  satisfy  the 

hypothesis  of  the  theorem.  Therefore,  by  blending  Mackey's  and  Stone's  theorem, 

the  following  commutative  diagram  results: 

(h,  e (u) . u (£) ) 


Figure  1 


Remarkably,  the  diagram  reveals  that  an  arbitrary  URS  is  equivalent,  under  a 
memoryless,  time  invariant,  unitary  transformation— a uniform  resolution  space 
insomorphism,  to  an  L2  space.  Distinctions  among  the  spaces  for  a fixed  G, 
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therefore,  depend  only  on  the  cardinality  of  the  space,  K.  Mixing  this 

equivalence  with  a recent  result  of  Masanl  and  Rosenberg10,  gives  the  desired 

structure— i.e.,  “time  invariance"  maps  to  multiplication. 

APPLICATION  OF  THE  MASANI -ROSENBERG  RESULT 

This  section  begins  with  the  result  of  the  above  mentioned  authors.  The 

theorem  is  not  stated  In  its  general  form10,  but  is  restricted  to  a group  J,  a 

Hilbert  space,  K,  and  the  Haar  measure,  m. 

Let  T be  a closed,  single-valued,  linear  operator  with  dense  domain  on 

L2(J,K,m),  such  that  T commutes  with  the  operation  of  multiplication  by  the 

characteristic  function  of  a Borel  set— i.e., 

X-F  e Tx«.  For  all  B In  £ 

^ ^ (th*  a-*lg»bra  of  all  Sorel 

secs  of  J). 

The,  there  exists  a measurable  function,  T on  J,  whose  values  are  operators 
on  K,  such  that 

(Tf)(j)  • T(J}f(J)  . j In  j. 

This  tneorem  applies  to  function  space.  Thus,  to  verify  the  sought  after 
properties  on  the  abstract  URS,  we  fist  apply  the  Mackey  Transforms  as  In  Figure  1, 
redrawn  below  for  simplicity. 

(H.  £ (u),  u (£)) 

(L^CG.X),  Xg * ° )(“  ' ’ Xg)> 

Our  discussion  dwells  upon  two  types  of  operators  in  the  abstract  URS,  memory- 
less operators  and  time  invariant.  First,  we  consider  the  time  Invariant  case. 

Let  T be  a closed,  linear,  single-valued  time  invariant  operator  on 
(H,E(U),  U(E)) . Recall  that  T is  time  Invariant  if  E(B)T  ^TE(B).  Under  the 
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Mackey  Transform,  M (which  we  term  the  Mackey  frequency  transform),  we  have  an 

a 

equivlent  statement  In  (L2(G,K),  oY,  X0),  as  follows: 

*btm  £ 

* 

where  T£  Is  the  Image  of  T under  the  M transformation.  • Clearly,  the  conditions 
of  the  Masanl -Rosenberg  theorem  are  satisfied.  Thus,  there  exists  a mapping,. 

A A 

T:&+K,  such  that 

(Th)(y)  - f(Y)h(y) 

A A 

for  all  y In  G and  h In  L2(G,k).  This  says  that  time  Invariant,  closed,  linear, 
single-valued  operators  on  an  abstract  URS  are,  as  hoped,  multiplications  In 
the  "frequency  doma1n"--1 .e. , In  (L2,(G,K) ,oY,  Xg). 

Now,  let  T be  a memoryless,  linear,  closed,  single-valued  operator  on 
(H,E,U) . Recall  that  T Is  memoryless  If  E*T  ■ TEt.  Thus,  by  reasoning  similar 
to  the  time  Invariant  case,  the  Image  of  T under  the  Mackey-time  transform,  M, 
commutes  with  Xg  In  (L2(G,K),  Xg,  c*).  Thus,  It  Is  equivalent  to  a multipli- 
cation by  the  Masanl -Rosenberg  theorem. 

This  structure,  then,  shows  that  certain  operators  on  an  abstract  URS  have 
the  "right"  properties.  It  Is  Interesting  to  note  that  there  Is  a duality  In- 
herent In  this  formulation.  The  presence  of  a Mackey  "time- transform"  and  a 
corresponding  "frequency-transform"  Is  apparently  necessary  for  the  cohesive- 
ness of  the  theory. 

CONCLUSIONS 

The  above  Ideas  assimilate  past  theories  In  a number  of  ways.  The  theory 

generalizes  the  Falb-Freedman-Anton  work  because  of  the  abstract  setting  and 

because  It  Is  valid  for  a larger  class  of  operators.  Clearly,  there  is  no  re- 

strlction  scalar-valued  frequency  responses  as  in  (9)  and  (13).  Moreover,  It 

21 

circumvents  the  existence  questions  associated  with  Saeks'  work.  In  fact,  as 


t 
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In  (21),  given  appropriate  conditions,  a multiplication  on  a function  space  can 
be  viewed  as  an  Integral  .over  the  spectral  measure,  defined  via  multiplication 
by  X,  l.e. 

T - / T(«s)<ixU). 

Hence,  the  pre-image  of  T under  the  Mackey  frequency  transform  In  (H,E,U)  takes 
the  form 

T - / T(r)«il(r). 
as  was  specified  In  (21). 

In  a private  conversation  with  one  of  the  authors,  Desoer  raised  a question 
about  the  fact  that  differential  operators  satisfied  the  definition  of  memory- 
lessness  in  an  abstract  setting.  The  question  caused  some  doubt  In  our  minds  as 
to  the  appropriateness  of  the  definition.  This  note  gives  a parltal  answer.  In 
that  closed,  memoryless  operators  are  multiplications.  Hence,  the  apparent 
pathology  noted  by  Desoer  can  only  arise  in  the  case  of  nonclosed  operators. 
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CHAPTER  I 
INTRODUCTION 

One  technique  which  has  been  developed  in  recent 
years  for  the  operator  theoretic  study  of  systems  is  the 
use  of  resolution  space  [1] . The  basic  motivation  behind 
the  development  of  resolution  space  was  to  overcome  the 
impossibility  of  defining  time-based  concepts  such  as 
causality  in  the  Hilbert  and  Banach  spaces  which  are  the 
setting  for  classical  operator  theory. 

Resolution  space  techniques  have  been  very  successful 
in  achieving  the  goal  of  including  time-based  concepts  in 
operator  theoretic  systems  theory,  but  the  concept  of  time 
which  has  been  used  is  classical  in  nature.  It  has  been 
known  since  the  early  part  of  this  century  that  an  accurate 
model  of  the  physical  universe  must  be  based  on  the  concept 

9 

that  space  and  time  are  intimately  connected.  This  is  the 
central  thesis  of  the  theory  of  relativity  formulated  by 
Einstein  [2] . Classical  resolution  space  techniques  ignore 
this  connection  between  space  and  time,  and  thus,  it  might 
be  suspected  that  an  extension  of  the  resolution  space 
concept  to  include  the  constraints  posed  by  the  theory  of 
relativity  could  offer  new  insights. 

Due  to  the  difficulty  of  merging  the  operator  theory 
required  by  resolution  space  and  the  theory  of  differenti- 
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able  manifolds  required  by  general  relativity,  the  relativ- 
istic resolution  space  theory  is  developed  only  for  the  case 
comparable  to  the  classical  resolution  space  development 
based  on  Hilbert  spaces  such  as  H - l2(G,k,u),  the  Hilbert 
space  of  functions  defined  on  an  ordered,  locally  compact, 
abelian  group  G which  take  values  in  a Hilbert  space  K,  and 
which  are  square  integrable  relative  to  a Borel  measure  y. 

In  H,  we  can  define  a spectral  measure  E by  multiplication 
by  the  characteristic  function 

1.1.  [E (A) f 1 (s)  - XA(s)f (s) 

for  each  Borel  set  A.  Given  the  spectral  measure  E,  we  can 
define  a resolution  of  the  identity  by 

1.2.  Efc  - E f°°, t ) , t e G. 

In  this  case,  Et  reduces  to  a family  of  truncation  operators 

f (s)  ; s < t 

1.3.  (Etf) (s)  - 

0 ; s > t . 

We  also  have  a resolution  of  the  identity  Et  defined  by 

1.4.  Efc  • I - Et. 

The  pair  (H,E)  is  the  classical  L2  resolution  space  in  which 
the  resolution  of  the  identity  Efc  allows  the  introduction  of 


the  desired  time-based  concepts. 


We  define  space-time  to  be  a pair  (M,g)  where  M is  a 
real,  four-dimensional,  connected  C*  Hausdorff  manifold,  and 

t 

g is  a globally  defined  C*  tensor  field  of  type  (0,2)  which 
is  nondegenerate  and  Lorentzian.  Then  analogously  to  the 
above  definition  for  classical  L2  resolution  space,  we  can 
define  relativistic  L2  resolution  space  to  be  the  pair  (H,E) 
where  H - L-,(M,K,y)  is  the  Hilbert  space  of  functions  de- 
fined on  M with  values  in  a Hilbert  space  K,  and  which  are 
square  integrable  with  respect  to  a Borel  measure  y,  and  E 
is  a spectral  measure  defined  by  multiplication  by  the 
characteristic  function 

1.5.  [E (A) f ) (s)  - XA(s)f (s) 

for  each  Borel  set  A. 

In  order  to  proceed  further,  it  is  necessary  to  find 
some  counterpart  to  the  resolution  of  the  identity  induced 
in  the  classical  case  by  the  spectral  measure  E.  We  don't 
have  a resolution  of  the  identity  in  the  relativistic  case 
since  the  manifold  M isn't  ordered.  It  is  still  possible  to 
parallel  the  classical  resolution  space  development  if  we 
first  define  the  past  and  future  of  any  point  x in  M. 


Although  some  care  must  be  taken  to  obtain  a precise  defi- 
nition, the  past  of  a point  x is  essentially  the  set  of  all 
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points  in  M which  could  have  sent  a signal  in  the  past 
traveling  at  a speed  less  than  or  equal  to  the  speed  of 
light  which  could  be  received  by  x in  the  present.  Simi- 
larly, the  future  of  x is  the  set  of  all  points  in  M which 
could  receive  a signal  sent  by  x traveling  at  a speed  less 
than  the  speed  of  light.  Then  the  family  of  projections  E5 
is  defined  by 


1.6. 


E[J  (x)] 


where  j”(x)  is  the  past  of  x,  and  similarly,  the  family  of 
projections  E^  is  defined  by 


1.7. 


Ex  - E[,T  (x)  ] 


where  J (x)  is  the  future  of  x. 

In  addition  to  the  fact  that  Ex  and  Ex  don't  form 
resolutions  of  the  identity,  we  also  have 


1.8. 


Ex  + Ex  ? I (Since  J+(x)U  j“(x)  t M)  , 


whereas  in  the  classical  case 


1.9. 


Et  + Et  - 1 


The  lack  of  an  order  on  M,  and  the  noncomplementary  nature 
of  Ex  and  Ex  combine  to  make  the  relativistic  definitions  of 
concepts  such  as  causality  and  strict  causality  more  com- 
plicated, and  definitions  which  were  equivalent  in  the 
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classical  case  are  no  longer  equivalent  in  the  relativistic 
case.  It  also  turns  out  that  strictly  causal,  strictly 
anticausal,  and  memory less  operators  are  no  longer  enough 
for  decomposition  of  an  arbitrary  operator,  and  the  totally 
new  concept  of  a spacelike;  operator  must  be  introduced.  The 
increased  complexity  of  the  relativistic  case  also  prevents 
several  of  the  classical  theorems  from  carrying  over  into 
the  relativistic  setting. 

The  most  familiar  space-time  is  the  Minkowski  space -time 
of  special  relativity  ( [2]  , [3]  , [4] ) . Minkowski  space-time 

4 

is  the  manifold  R with  a flat  Lorentz  metric  g.  If  x = 

(x1,x2,x3,x4)  and  y = (y1,y2*y3y4)  are  two  points  in  Min- 
kowski space-time,  then 

1.10.  g(x,g)  * + x2y2  + x3y3  - x4y4. 

A nonzero  point  x is  said  to  be  timelike  of  g(x,x)<0, 
spacelike  if  g(x,x)>0,  and  null  if  g(x,x)  = 0.  For  a given 
point  x,  all  the  points  separated  from  x by  a timelike  or 
null  distance  form  a hypercone  called  a lightcone.  If  we 
assume  that  all  timelike  distance  vectors  can  be  classified 
as  either  future-directed  or  past-directed,  then  the  light- 
cone  can  be  divided  into  two  parts.  Points  separted  from 
x by  a future-directed  timelike  distance  are  said  to  be  in 
the  future  of  x and  points  separated  from  x by  a past- 
directed  nonspacelike  distance  are  said  to  be  in  the  past  of 
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x.  All  other  points  are  said  to  lie  in  the  spacelike  region 
about  x.  With  two  spatial  dimensions  supressed,  the  light- 
cone  of  a point  x is  shown  in  Figure  1. 

In  the  next  chapter,  we  develop  a relativistic  reso- 
lution space  theory  for  the  special  case  in  which  the  space- 
time  manifold  is  a two-dimensional  version  of  Minkowski 
space-time.  We  first  define  a special  set  of  lines  called 
null  lines.  These  lines  are  essentially  the  paths  along 
which  a light  ray  would  travel  in  our  space-time.  Then  the 
past  and  future  of  a null  line  are  defined,  and  these  con- 
cepts are  used  to  define  a special  class  of  sets  called 
diamond  sets.  The  diamond  sets  are  shown  to  form  a semi- 
ring and  this  permits  them  to  be  used  to  establish  an  in- 
tegration theory  paralleling  that  used  in  the  classical 
resolution  space  development  ( [1] , [5] ) . 

Next,  causal,  anticausal,  and  memoryless  operators  are 
defined,  and  then  shown  to  have  properties  similar  to  those 
in  the  classical  development  [1] . 

After  this,  strictly  causal  and  strictly  anticausal 
operators  are  defined,  and  once  again,  the  development  in 
the  classical  case  is  paralleled  [1] . 

Finally,  the  decomposition  theorems  for  an  arbitrary 
operator  are  stated  and  proved  [1].  In  order  to  obtain  a 
complete  decomposition,  a new  type  of  operator  called 
spacelike  is  introduced.  This  operator  essentially  takes 
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care  of  the  spacelike  region  ignored  by  the  causal  and 
anticausal  operators. 

The  conclusion  discusses  the  problems  involved  in 
extending  the  two-dimensional  theory  to  a more  general  case. 
A possible  technique  for  developing  a general  theory  is 
indicated,  and  we  state  a conjectured  property  of  the 
spacelike  operators. 
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CHAPTER  II 

THE  TWO  DIMENSIONAL  CASE 

Our  space-time  will  be  an  ordered  pair  (M,g)  where  M 
is  the  connected,  two-dimensional  (one  space  and  one  time 

OB  2 

coordinate) , Hausdorf f C manifold  R^  and  g is  the  usual 
flat  Lorentz  metric  used  in  special  relativity,  i.e.,  for 
( x^ , t ^ ) , (X2 , tj) eM , 

2.1.  g [ (x1,t1) , (x2,tfc) ] - *1*2  “ tit2 * 

2 

This  space  is  essentially  R with  an  "inner  product'*  defined 
by  g.  Following  Penrose  [6] , a non-zero  tangent  vector  X is 
said  to  be  time like  if  g(X,X)<0,  spacelike  if  g(X,X)>0,  and 
null  if  g(X,X)  » 0.  A C^*  curve  in  M is  called  time  like 
(spacelike  , null)  if  the  tangent  vector  to  the  curve  at 
each  point  is  timelike  (spacelike,  null) . A curve  will  be 
called  non-spacelike  if  the  tangent  vector  at  each  point  of 
the  curve  is  either  timelike  or  null. 

We  now  assume  that  we  can  divide  the  non-spacelike 
vectors  at  each  point  in  M into  two  groups  which  we  will 
call  future-  and  past-  directed  non-spacelike  vectors. 
Essentially,  a non-spacelike  vector  is  future-directed  if  it 
makes  an  angle  of  45°  to  135°  with  the  x-axis,  and  it  is 
past-directed  if  it  makes  an  angle  of  225°  to  315°  with  the 
x-axis. 
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Let  x ■ (x,t)  be  a point  of  M.  The  future  of  x will 
be  denoted  J+(x)  and  is  defined  to  be  the  set  of  all  points 
of  M which  can  be  reached  from  x by  a future-directed  time- 
like curve,  i.e.,  a curve  whose  tangent  vectors  are  all 
future-directed  timelike  vectors.  J+(x)  doesn't  include 
x.  The  past  of  x will  be  denoted  J~(x)  and  is  the  set  of 
all  points  of  M which  can  be  reached  from  x by  a past- 
directed  non-spacelike  curve,  i.e.,  a curve  in  M whose 
tangent  vectors  are  all  past-directed  non-spacelike  vec- 
tors. J"(x)  does  include  x.  For  an  illustration,  see 
Figure  2. 

It  is  easily  seen  that  the  boundaries  of  the  past  and 
future  of  a point  are  lines  with  slopes  of  + 1.  Since  these 
lines  are  so  important,  they  will  be  given  a name,  null 
lines.  The  future  of  a null  line  is  defined  to  be  the  set 
of  all  points  of  M which  lie  in  the  future  of  some  point  on 
the  null  line.  The  past  of  a null  line  is  defined  to  be  the 
set  of  all  points  of  M which  lie  in  the  past  of  some  point  on 
the  null  line.  See  Figure  3. 

In  order  to  carry  through  the  integration  theory  in  our 
space-time  setting,  we  need  to  define  the  class  of  sets 
which  will  be  used  to  partition  M.  A diamond  set  will  be 
defined  to  be  the  intersection  of  the  pasts  or  futures  of 
any  finite  set  of  null  lines,  the  empty  set,  or  the  whole 
space  M.  The  future  (past)  of  a diamond  set  is  the  union  of 
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the  futures  (pasts)  of  all  the  points  contained  in  the 
diamond  set  and  is  denoted  J+ (D) [J~ (D) ] . All  possible 
diamond  sets  are  shown  in  Figure  4. 

Our  first  theorem  concerns  the  suitability  of  the 
diamond  sets  for  performing  the  integrations. 

2.1  THEOREM  The  diamond  sets  form  a semiring. 

Proof:  To  establish  this,  we  need  to  show  that 

i)  the  intersection  of  two  diamond  sets  is  a diamond 
set,  and 

ii)  the  set  difference  of  two  diamond  sets  is  a dis- 
joint union  of  diamond  sets. 

i)  is  immediately  obvious  from  the  definition  of  a 
diamond  set.  ii)  is  apparent  from  the  diamond  sets  pictured 
in  Figure  4 and  the  fact  that  all  lines  bounding  the  diamond 
sets  have  a slope  of  + 1. 

In  order  to  keep  the  analogy  with  the  classical  case, 
we  would  like  to  be  able  to  write  the  past  and  future  of 
each  diamond  set  as  the  past  of  a single  point.  With  each 
null  line  L,  associate  two  additional  points  ur  and  lr  . 

will  be  called  the  upper  point  of  L and  has  the  property 
that  the  past  of  u_  is  the  past  of  L.  The  future  of  u_  will 
be  the  empty  set.  I L will  be  called  the  lower  point  of  L 
and  has  the  property  that  the  future  of  is  equal  to  the 
future  of  L.  The  past  of  is  considered  to  be  the  empty 
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set.  We  also  add  two  more  points,  and  -«.  The  past  of 
and  the  future  of  -<*»  are  M,  and  the  future  of  +«*>  and  the 
past  of  -*  are  the  empty  set.  These  added  points  are  purely 
a notational  convenience,  and  as  such,  they  have  no  relation 
to  the  points  of  M.  We  have  not  yet  been  able  to  discover 
whether  or  not  there  is  a topology  which  would  continuously 
extend  the  metric  to  include  these  points.  With  the  inclu- 
sion of  these  extra  points,  it  is  now  possible  to  write  the 
past  and  future  of  any  diamond  set  as  the  past  and  future  of 
two  unique  points.  For  an  example,  see  Figure  5. 

Now  let  L2(M)  be  the  Hilbert  space  of  L2  functions 
defined  on  M.  If  A is  a subset  of  M,  we  define  the  pro- 
jection 


2.2.  E (A) : L2(M)  * L2(M) 


by 

f (x) , xeA 

2.3.  E (A) f ( x)  - , f £ L,(M)  . 

0 , xfA  Z 


Analogously  to  the  classical  case,  we  define  the  projections 
EX:  L2  (M)  * L2(M)  and  Ex:  Lj(M)  -*■  L2(M)  for  x e M by 

2.4.  EX  - E [ j”  ( x) ] 

and 

Ex  - E[J+(X)). 


2.5. 


Let  T be  an  operator  on  LjW  . We  say  T is  causal  if  it 
satisfies  the  two  following  conditions: 

2.6.  EX  T ■ EX  T EX  V x e M 
and 

2.7.  T Ex  » Ex  T Ex  V x e M. 

Notice  that  in  the  classical  case,  these  two  causality 
conditions  imply  each  other.  This  is  a result  of  the  fact 
the  ranges  of  Ex  and  Ex  are  complementary  subspaces  in  the 
classical  case.  This  is  not  the  case  in  the  relativistic 
setting,  and  both  parts  of  the  definition  are  required  to 
give  a physically  desirable  interpretation  of  causality. 

We  say  that  T is  anticausal  if 

2.8.  E T - Ex  T Ex  V X 6 M 

2.9.  T EX  » EX  T EX  V X £ M. 

Let  D(u,A)  denote  a diamond  set  D with  J+(D)  = J+(A) 
and  J~(D)  ■ J-(u).  We  have  the  following  theorem. 

2.2  THEOREM  T:  Lj  (M)  -*■'  L2(M)  is  causal  if  and  only  if 

2.10.  E [D (u , 1) ] T - E [D (u, A) ] T EU 
and 


2.11.  T E [D  (u , A)  ] * E^T  E[D(u,A)] 


for  all  diamond  sets  D(u,Z). 

Proof:  Suppose  T is  causal.  Then  EUT  » EUTEU  and  TE^  • E ^TE ^ 
It  is  easily  seen  that  E[D(u,Z)]Eu  » [D(u,Z)]  and  EaE[D(u,Z)] 

■ E [D  (u , Z) ] . Then  E[D(u,Z)]EUT  » E [D (u, Z) ] EUTEa  or  E[D(u,Z)]T 
» E [D (u, Z) ] TEa.  Similarly  TE[D(u,Z)]  - E^TE [D (u, Z) ] . 

Now  suppose  TE [ D ( u , Z ) ] ■ E^TE  [D  (u, Z)  } and  E[D{u,Z)]T  =» 

E [D (u , Z) ] TEU  for  all  diamond  sets.  Let  x be  any  point  of  M. 
Then  J+(x)  and  J”(x)  are  diamond  sets.  It  is  also  easily 
seen  that  J+[J+(x)]  ■ J+(x)  and  j”[J~(x)]  * J~(x).  From 
this,  it  follows  that  EXT  * EXTEX  and  TEX  * E ^TE x , 

We  obtain  a similar  theorem  for  anticausal  operators. 

2.3  THEOREM:  T:  L2(M)  L2(M)  is  anticausal  if  and  only 

if 

2.12.  TE [D (u, Z) ] - EUTE[D(u,Z)]  , 
and 

2.13.  E[D(u,Z)]T  » E[D(u,Z)  ]TEr 

An  operator  T:  L2(M)  ■+  L2(M)  is  said  to  be  memory  less 
if  T is  both  causal  and  anticausal.  We  have  the  following 
theorem. 

2.4  THEOREM:  T is  memoryless  if  and  only  if 

2.14.  E[D(U,Z)]T  * TE[D(u,Z)] 
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for  every  diamond  set  D(u,l). 

Proof:  First  assume  T is  memoryless.  Then  T is  both  causal 
and  anticausal.  T causal  implies  that 

E[D(u,£)]T  - E[D(u,£) ]TEU  , and 

TE  [D  (u,  £)  ] - E^TEIDU,*)  ] , 

and  T anticausal  implies  that 

E[D(u,£) ]T  * E[D(u,£) JTE^  , and 

TE [D ( u , £)]  = EUTE[D(u,£) ] . 

Then  E[D(u,i)]T  = [E (D (u, £) ) EU  * [E (D (a, £ ) ) TE^] EU 

- E ( D ( u , £ ) ) TE ( D ( u , £ ) ) , 


and 

TE [ D ( u , £ ) ] * E^ [TE (D  (u, £) ) ] 

* E4IEuTE{D(u/£)  ) ] 

* E(D(u,£))TE(D(u,£) ) . 

Hence  E[D(u,£)]T  = TE[D(u,£)]  for  all  diamond  sets  D(u,£). 

Now  suppose  E[D(u,£)]T  * TE[D(u,£)]  for  all  diamond 
sets  D (u, £) . 

Then  E [D (u, £) ] T » E [D (u , £) ] TE [D (u , £) ] 


t 
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and  E[D(u,f)  ] TEU  - E [ D ( u , l)  ] TE  [ D { u , & ) ]Ea 

- E [D (a, i) ] TE [D (a, A) ] - E[D(u,i)]T. 

Similarly  E^TElDtu,!)]  ■ TE[D(u,i)].  Hence  T is  causal. 
Similarly  T is  anticausal.  Therefore  T is  memory less. 

The  classes  of  operators  which  have  been  defined  are  as 
well-behaved  as  in  the  classical  case. 

2.5  THEOREM;  The  set  of  causal  (anticausal,  memoryless) 
bounded  linear  operators  from  a Banach  algebra  with  identity 
which  is  closed  in  the  strong  operator  topology  of  the 
algebra  of  all  linear  bounded  operators. 

Proof:  The  proof  will  be  presented  for  the  causal  case  only 

since  the  proof  for  the  anticausal  case  is  very  similar,  and 
the  result  for  the  memoryless  case  follows  from  the  fact 
that  the  intersection  of  two  Banach  algebras  is  also  a 
Banach  algebra. 

If  T and  S are  causal,  then 

EXTS  - (EXTEX)S  = EXT(EXS)  = EXT(EXSEX) 

» (EXTEX)SEX  * EXTSEX  . 

Similarly,  TSEX  » ExTSEx,  and 

hence  ST  is  causal.  If  we  take  the  sum  of  S and  T,  we  have 
EX(S+T)«  EXS  + EXT  - EXSEX  + EXTEX  = EX(S  + T)EX,  and 
(S  + T)EX  - Ex(S  + T)EX. 

Thus  S + T is  causal.  Also,  EXI  * EXEXI  = EXIEX  and  IEx 
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» e^IE^ . This  proves  the  identity  is  causal.  Finally,  if 

is  a sequence  of  causal  operators  converging  strongly  to  T, 

i.e.,  lim  T.f  * Tf  for  all  f,  then  since  Ex  is  bounded, 
i-~»  1 

EXTf  = EX [lim  T f]  » lim  EXT . f 
±-+<0 

* lim  EXT.EXf  » EX [lim  T . EXf ] = EXTEXf . 

Similarly,  TExf  = E^TE^f  for  all  f.  Hence  T is  causal. 

We  have  the  following  theorem  for  the  adjoint  of  a 
causal  operator. 

2.6  THEOREM:  An  operator  T is  causal  if  and  only  if  T*  is 
anticausal . 

Proof:  Suppose  T is  causal.  Then  EXT  * EXTEX.  Taking  the 

adjoint  of  both  sides,  we  have  (EXT) * = (EXTEX) * which 
reduces  to  T*EX  =»  EXT*EX.  Similarly,  E^T*  = ExT*Ex«  Hence 
T*  is  anticausal.  The  converse  is  similar. 

Since  a memoryless  operator,  is  both  causal  and  anti- 
causal, then  so  is  its  adjoint  and  hence  the  adjoint  of  a 
memoryless  operator  is  also  memoryless. 

At  the  moment,  nothing  can  be  said  about  the  inverse  of 
a causal  operator  in  the  relativistic  case.  The  results 
from  the  classical  case  don't  carry  over,  partly  because  the 
ranges  of  the  projections  aren't  complementary  subspaces, 
and  partly  because  the  condition  Exf  * Exg  =>  ExTf  * 
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ExTg,  v x e M,  f,geL2(M),  is  not  enough  to  insure  causality 
for  T in  the  relativistic  case. 

We  cone  now  to  the  extension  of  the  integrals  of 
triangular  truncation  to  the  relativistic  setting.  It  was 
for  this  purpose  that  the  diamond  sets  were  introduced. 

Since  they  form  a semiring,  they  can  be  used  to  partition 
M so  that  the  integrals  can  be  defined  [5].  The  upper 
Cauchy  integral  of  an  operator  valued  function  f on  M is 
defined  by 

, niP> 

2.15.  UC  / f(x)dE(x)  - lim  £ f (u . ) E [D . (u . , Z . ) ) 

peP  i«l  1 111 

where  the  limit  exists  in  the  uniform  topology  over  the  net 
of  all  partitions  of  M into  diamond  sets  Simi- 

larly, the  lower  Cauchy  integral  is  defined  by 

2.16.  LC  / f ( x)  dE  ( x)  - lim  f ( Z.  ) E [D.  (a  .,£.)]  . 

peP  i-1  1 ill 

These  integrals  can  also  be  defined  with  the  measure  on  the 
left,  or  over  a portion  of  M instead  of  all  M.  We  can  also 
define  the  strong  Cauchy  integrals  SUC / and  SLC / by  taking 
the  limit  in  the  strong  operator  topology.  We  then  have  the 
following  theorem  relating  causality  to  these  Cauchy  inte- 
grals . 


C 


2.7  THEOREM:  The  following  are  equivalent  for  a linear 
bounded  operator  T. 
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i)  T is  causal. 

2.17.  ii)  UC  / dE(x)TEX  - LC  / E^TdE ( X)  - T. 

2.18.  iii)  SUC  / dE(x)TEX  - SLC  / E xTdE ( X ) « T. 

Proof:  i)  ->  ii) . 

If  T is  causal,  then  E[D(u,£)]TEX  » E[D(u,£)]T  and  ExTE[D(u,l)] 
- TE  [D (u,  £.)  ] for  all  xeM.  Hence  for  any  partition  of  M into 
diamond  sets, 

n(p) 

T = IT  = [lim  l E [D . (u .,£.)]] T 
peP  i-1  111 

n(p) 

* lim  I E[D. (u. ,4. ) ]T 
peP  i-1  iii 

n(p)  u. 

- lim  l E [D . (u. ,1.) ]TE 

peP  i-1  i 1 i 

- UC  / dE(x)TEX. 

Similarly,  T - LC  / E^TdE (x) . 

ii)  ->  iii) . 

This  follows  immediately  from  the  fact  that  uniform 
con  .egence  implies  strong  convergence. 

iii)  ->  i. 

We  have  SUC  / dE(x)TEX  - SLC  / E ^TdE ( x ) - T.  Then  E^T  » 

y 

E SUC/ydE  ( x)  TEX  - SUC/j-(t/)dE(x)TEX  - SUC/j-  {(j)  dE  (X)  TeV' 
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- [SUC/T-,  ,vdE(x)TEX]E{/.»  E^TE^.  Similarly,  TE  , ■ E ,TE,. 

j y y y 

Hence  T is  causal. 

A similar  theorem  is  also  true  for  anticausal  operators. 

2.8  THOEREM:  The  following  are  equivalent  for  a bounded 
linear  operator  T. 

i)  T is  anticausal. 

2.19.  ii)  UC/E*TdE(x)  - LC/dE(x)TEx  - T. 

2.20.  iii)  SUC/E*TdE(x)  « SLC/dE(x)TEx  - T. 

For  an  integral  representation  of  memoryless  operators, 
we  need  to  define  the  diagonal  Cauchy  integral 
, niP> 

2.21.  C / dE ( x ) TdE ( X ) -lim  £ E [D. (u. , l . ) ]TE [D, (u . , l . ) ] 

peP  i-1  111  ill 

where  the  limit  is  taken  in  the  uniform  topology  over  the 
net  of  all  partitions  of  M into  diamond  sets.  We  can  also 
define  the  strong  diagonal  Cauchy  integral  SC/dE(x)TdE(x) 
by  taking  the  limit  in  the  strong  topology.  We  then  have 
the  following  theorem  for  memoryless  operators. 

2.9  THEOREM : For  a bounded  linear  operator  T,  the  fol- 
lowing are  equivalent. 

i)  T is  memoryless. 

2.22.  ii)  C /dE ( x ) TdE ( x ) - T. 
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iii)  SC /dE ( X ) TdE ( X ) - T. 


Proof:  Similar  to  the  proof  of  the  theorem  for  causal 
operators. 

We  now  define  the  notion  of  strict  causality.  A 
bounded  linear  operator  T is  said  to  be  strictly  causal  if 

2.24.  LC/dE (x) TEX  « UC/ExTdE(x)  » T. 

A bounded  linear  operator  is  said  to  be  strongly  strictly 
causal  if 

2.25.  SLC/dE(x)TEX  » SUC/ExTdE(x)  = T. 

The  strictly  causal  case  is  different  from  the  causal  case 
in  that  strict  causality  implies  strong  strict  causality 
but  not  conversely  (see  [1]). 

In  order  to  characterize  the  relationship  between  the 
strictly  causal  and  causal  operators,  we  need  to  define  the 
following  integrals. 

v n(p)  Ui  *1 

2.26.  SfdE(x)TEX  - lim  ■ > E[D. (u. , £ . ) ]T(E  1 ” E ; . 

peP  i-1  1 1 1 

n(D) 

2.27.  S/EvTdE<X)  - lim  l (E9  - E(J  ) TE  [D . (u.  , l . ) ] . 

x oeP  i-1  x i ui  iii 

The  limit  is  taken  in  the  uniform  topology  over  the  net  of 
all  partitions  of  M into  diamond  sets.  We  can  also  define 
the  integrals 


SSjdEix)TEX  and  SS/ExTdE(x) 


where  the  1/  !t  is  taken  in  the  strong  topology.  We  then 
have  the  following  theorem  for  strictly  causal  operators. 


2.10  THEOREM:  A bounded  linear  operator  T is  strictly 
causal  if  and  only  if  T is  causal  and 

2.28.  S/dE(x)TEX  - ff/f  dE(x)  . 

Proof:  First  suppose  that  T is  strictly  causal. 

EXT  * E^C/ydEfi/JTE^  - LC/ j”  ( x)  dE  ( tj)  TE^ 

- LC/j-  ix)dEl{/)TE^EX  • EXTEX. 

Similarly  TEX  « E^TE^ . Hence  T is  causal.  T being  causal 
implies  that 


T » UC/dE ( x) TEX 

n(p)  u. 

- lim  I E [D . (u. ,1.) ]TE 

peP  i-1  111 

n (p)  u.  i 

- lim  V E [D . (u . , & . ) ] T(E  1 + [E  1 “ E X] ) 
peP  i-1  ill 

n (p)  i 4 n(p) 

- lim  ) E [D  . (u  . , 4 . ) ] TE  1 + lim  I E[D.  (u.  ,4.  ) ]' 

peP  i-1  i i i peP  i-1  ill 

u.  4 . 

[E  -E  ] 

« LC/dE(x)TEX  ♦ s/dE(x)TEX 
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Hence  s/dE(x)TEx  ■ 0.  Similarly  S/ExTdE(x)  ■ 0. 

By  reversing  the  above  argument,  it  can  be  seen  that  T 
causal  and  s/dE(x)TEx  » s/ExTdE(x)  - 0 implies  that  T is 
strictly  causal. 

We  have  a similar  theorem  for  strongly  strictly  causal 
operators. 

2.11  THEOREM : A bounded  linear  operator  T is  strongly 
strictly  causal  if  and  only  if  T is  causal  and 

2.29.  SS/dE(x)TEX  » SS/ExTdE(x)  =»  0. 

There  are  also  similar  theorems  for  strictly  anticausal 
and  strongly  strictly  anticausal  operators.  We  first  need 
to  define  the  following  integrals. 

n (p)  u.  1. 

2.30.  S / E xTdE  ( x ) « lim  £ [E  -E  j TE  [D . (u.  , £.)  ] . 

peP  i»l 

n^p) 

2.31.  S /dE  ( x ) TE  = lim  l E[D.  (u.  , l .)  ]T[E„  -E((  ]. 

“X  peP  i-1  111  i Ui 

Again,  the  limit  is  taken  in  the  uniform  topology  over  the 
net  of  all  partitions  of  M into  diamond  sets.  Then  we  have 
the  following  theorem  which  is  proved  in  the  same  manner  as 
the  corresponding  theorem  on  strictly  causal  operators. 

2.12  THEOREM : A bounded  linear  operator  T is  strictly 
anticausal  if  and  only  if  T is  anticausal  and 


* 
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2.32.  S/E*TdE(x)  » S/dE(X)TEx  - 0. 

If  we  take  limits  in  the  strong  topology,  we  obtain  the 
following  theorem  for  strongly  strictly  anticausal  operators. 

2.13  THEOREM?  A bounded  linear  operator  T is  strongly 
strictly  anticausal  if  and  only  if  T is  anticausal  and 

2.33.  SS/EXTdE(x)  * S5'/dE(x)TEx  - 0. 

We  have  the  following  theorem  for  the  space  of  strictly 
causal  (strictly  anticausal,  strongly  strictly  causal, 
strongly  strictly  anticausal)  operators. 

2.14  THEOREM : The  space  of  strictly  causal  (strictly 
anticausal,  strongly  strictly  causal,  strongly  strictly 
anticausal)  operators  forms  a Banach  space  which  is  closed 
in  the  uniform  operator  topology  of  the  space  of  all  bounded 
linear  operators. 

Proof:  The  proofs  in  all  four  cases  are  similar,  so  only 
the  strictly  causal  case  will  be  presented. 

These  operators  form  a Banach  space  since  they  are 
defined  by  a linear  equation  T » LC/dE(x)TEx.  Now  suppose 
T^T  where  the  T^  are  strictly  causal.  Since  the  T^  are 
strictly  causal,  they  are  causal,  and  hence  T is  causal.  We 
would  now  like  to  show  that  s/dE(x)TEx  » 0. 

For  any  e > 0,  choose  j such  that  | | T ^ - T| | < e/2  and  a 
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partition  of  M into  diamond  sets  D^(u^, 1^) , i - l,...,n  such 
that 


n 


u. 


n 


| | I'E[D. (u.,£.)]T. (E  X) | | < e/2 
i-1  1 1 1 3 

n u.  A. 

Then  II  [ E[Di(ui,£i)]T(E  -E  x)  I I « | | IeID.^,^) ] (T-T^ 

u.  if*1  n u.  £.  lml 

(E  1-E  X)  + lE[D. (u. ,£.)]T, (E  -E  X)  | | < sup | | E [D. (u . , £.) ] 

i»l  i* 

u.  j u.  £ . 

(T-T.ME  1-EXi)  | | + Syp| |E[D. (u. ,£.) ]T. (E  -E  l)  | | < e/2  + 

J i 1 J 

e/2  » e.  So  the  partial  sums  for  sfdE(x)TEx  converge  to 
zero.  Similarly,  s/ExTdE(x)  * 0.  Hence  T is  strictly 
causal . 

The  following  theorem  relates  the  strictly  causal  and 
the  strictly  anticausal  operators  to  each  other. 

2.15  THEOREM:  An  operator  T is  (strongly)  strictly  causal 

if  and  only  if  T*  is  (strongly)  strictly  anticausal. 

* n (p)  £.= 

Proof:  LC/dE ( x) TE  * lim  £ E [D. (u. , £ . )TE  . If  we  take 
n (p)  peP  i-1  H.1 

the  adjoint  of  T E [D . (u. , £ . ) ]TE  A,  then  we  have 
i-1  111 

n(p)  £ • 

v e XT*E [D. (u .,£.)] . Since  the  adjoint  is  a linear 

. ill 

i-l 

isometry  in  the  space  of  bounded  linear  operators  on  a 
Hilbert  space,  the  first  integral  converges  to  T if  and  only 
if  the  second  integral  converges  to  T*.  Using  a similar 
procedure  for  UC/Ex  TdE(x),  we  see  that  T is  strictly  causal 
if  and  only  if  T*  is  strictly  anticausal. 

In  order  to  state  and  prove  the  additive  decomposition 
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f 


* 


C 
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theorem  for  arbitrary  operators,  we  first  need  to  define  a 
new  class  of  operators  which  will  be  called  spacelike. 

For  a diamond  set  D(u,i)  we  define  the  projection 

2.34.  E*  - I-E*-Eu. 

With  this  definition,  we  say  that  an  operator  T is  spacelike 
if 

2.35.  TE  [D  (u , 2.)  ] - E*TE  [D (u,  l)  ] 
and 

2.36.  E[D(u,l)]T  - E [D (u, i) ] T£* 

for  all  diamond  sets  D(u,£). 

We  can  obtain  an  integral  characterization  of  spacelike 
operators  by  defining  the  following  integrals. 

, „ n(p) 

2.37.  sfE*TdE(x)  - lim  T E TE [D. (a . , i . ) ] . 

x peP  i-1  ui  111 

, x n(p)  l. 

2.38.  SfdE(x)TEx  « lim  l E[D. (u. , l . ) ]TE,  . 

x peP  i-1  111 

The  limit  is  taken  in  the  uniform  topology  over  the  set  of 
all  partitions  of  M into  diamond  sets.  We  can  also  take 
limits  in  the  strong  topology  in  which  case  we  obtain  the 
integrals  SSj ExTdE(x)  and  SS/dE{x)TEx.  We  then  have  the 
following  theorem. 
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2.16  THEOREM:  For  a bounded  linear  operator  T,  the  fol- 
lowing are  equivalent. 

i)  T is  spacelike. 

2.39.  ii)  s/E*TdE(x)  • S/dE(x)  TE*  - T . 

2.40.  iii)  SS/E*TdE(x)  - SS/dE(x)TE*  » T. 

Proof:  Similar  to  the  proof  for  the  causal  case. 

We  also  obtain  the  following  theorem  for  the  relation- 
ship between  a spacelike  operator  T and  its  adjoint  T*. 

2.17  THEOREM : T is  spacelike  if  and  only  if  T*  is  space- 
like. 

Proof:  TE[D(u,.l)]  - E*TE[D(u,i)  ] <=>E [D (a, fc)  ] T* 

■ E [D (u , l) ]T*eJ  and, 

E [D  (u , £) ] T - E[D(u, l) ]TE*<->T*E[D(u,4) ] 

« E^T*E [D (u , & ) ] . 

Finally,  we  obtain  the  additive  decomposition  theorem 
for  an  arbitrary  bounded  linear  operator. 

2.19  THEOREM:  Let  T be  an  arbitrary  bounded  linear  oper- 
ator. Then  T can  be  decomposed  as  T * C + A + S where  C is 


strictly  causal,  A is  strictly  anticausal,  and  S is  spaca 
like  if  and  only  if  the  integrals 
UC/ExTdE(x)  , LC/dE(x)TEX  , 

LC/EXTdE(x)  , UC/dE(x)TEx  , 

S/E^TdEU)  and  S/dE(x)TEx 

exist  and  satisfy 

2.41.  UC/ExTdE(x)  - LC/dE(x)TEX, 

2.42.  LC/EXTdE(x)  - UC/dE(x)TEx, 
and 

2.43.  S/EXTdE(x)  - S/dE(x)TEX  , 

in  which  case  the  decomposition  is  given  by 

2.44.  C - UC/ExTdE(x)  ■ LC/dE(x)TEX, 

2.45.  A - LC/EXTdE(x)  ■ UC/dE(x)TEx, 
and 

2.46.  S - S/EXTdE(x)  - S/dE(x)TEX. 

Proof:  First  assume  that  all  the  integrals  exist  and  sat 
isfy  the  required  relationships.  Then 
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n(p) 

T - UCjTdE(x)  - lira  l TE  [Di (ul^  l L)  ] 

peP  i»l 

n (p)  4.  4. 

- lira  I [E  1 + Eu  + Eu1]TElDi(ai4.)] 

peP  i-1  i i 

n(p)  4.  n(p) 

- lim  V E XTE  [D . {a,  , 4 . ) ] +lim  V TE [D. (u., 4.) ] 

peP  i-1  1 1 1 pePi-1  1 1 1 

n(p)  4. 

-►  lim  T E,,TE[D.  (a.  ,4.)  ] 
peP  i»l  i 

- LC/EXTdE(x)+UC/ExTdE(x)+S/E^TdE(x) . 

Similarly, 

T - LC/dE(x)T 

* UC/dE(x)TEx+LC/dE(x)TEX+S/dE(x)TEX 

We  have 

UC/E^CdEU)  * UC/  [UC/E^TdE  ( y)  1 dE  ( x) 

- UC/EzrdE(r)  * C. 

Similarly,  LC/dE(x)CEx  - C. 

Hence  C is  strictly  causal. 

Similarly,  A is  strictly  anticausal,  S is  space- 
like, and  therefore  this  is  the  desired  decomposition. 

Now,  suppose  that  the  desired  decomposition  exists. 

Then 


I 
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UC/ExTdE(x)  - UC/Ex(C+A+S)dE(x) 

- UC/ExCdEU)+UC/ExAdE(x)+UC/ExSdEU)  . 

n (p) 

UC/E  AdE(x)  - lim  l E AE[D. (a. ,1.)] 
x pcPi-1  ui  111 

n(p)  a. 

- lim  T E E XAE  [D . (u . , i. ) ] • 0. 

peP  i-1  ui  ill 

n (p) 

UC /E  SdE  ( x)  - lim  £ E„  SE  [D . (a  . . I . ) ] 
x peP  i-1  ui  ill 

n(p)  l. 

- lim  l E E,/SE(D.  (u.  ,4.)  ] - 0. 
p £P  i-1  ui  ai  ill 

Hence  UC/ExTdE(x)  - UC/ExCdE(x)  - C. 

Similarly  LC/dE(x)TEx  - C.  Thus  C is  of  the  required 
form.  In  the  same  manner,  A and  S can  be  shown  to  have 
the  required  form. 

We  have  a similar  theorem  for  additive  decomposition 
into  operators  defined  by  strongly  convergent  integrals. 

2.20  THEOREM:  Let  T be  an  arbitrary  linear  bounded 
operator.  Then  T can  be  decomposed  asT-C+A+S 
where  C is  strongly  strictly  causal,  A is  strongly 
strictly  anticausal,  and  S is  spacelike  if  and  only  if 
the  integrals 
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SUC/ExTdE(x)  , SLC/dE ( x) TEX , 

SLC/EXTdE(x)  , SUC/dE(x)TEx, 

SS/EXTdE(x)  and  SS/dE(x)TEX 
X ^ 

exist  and  satisfy 

2.47.  SUC/ExTdE(x)  = SLC/dE (x)TEX  , 

2.48.  SLC/EXTdE(x)  = SUC/dE(x)TEx  , 
and 

2.49.  SS/EXTdE(x)  = SS/dE(x)TEX, 

in  which  case  the  decomposition  is  given  by 

2.50.  C * SUC/ExTdE(x)  = SLC/dE ( x) TEX , 

2.51.  A “ SLC/EXTdE(x)  - SUC/dE ( x) TEX , 
and 

2.52.  S » SS/EXTdE(x)  - SS/dE(x)TEX  . 
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CHAPTER  III 
CONCLUSIONS 

Following  the  development  of  resolution  space  for  the 
two-dimensional  special  relativistic  case,  it  would  be 
nice  to  extend  the  development  to  the  four-dimensional 
special  relativistic  case,  and  then  to  the  general  rela- 
tivistic case.  However,  no  way  has  been  found  to  extend 
the  results  in  this  paper  to  even  the  three-dimensional 
special  relativistic  setting. 

One  major  problem  has  been  the  attempt  to  keep  a 
point  development  of  resolution  space.  The  lack  of  an 
ordering  on  the  space-time  manifolds  presented  difficulties 
even  in  the  two-dimensional  case.  By  using  diamond  sets 
it  was  possible  to  keep  a semblance  of  the  point  develop- 
ment, but  even  in  this  case  special  points  had  to  be 
introduced,  and  attempts  to  extend  the  space-time  topology 
to  include  the  extra  points  in  such  a manner  as  to  con- 
tinuously extend  the  metric  were  fruitless. 

We  attempted  to  extend  the  diamond  sets  to  higher 
dimensions,  but  we  were  unable  to  find  a higher  dimen- 
sional analogue  of  the  diamond  sets.  Hence,  it  would 
appear  that  the  diamond  sets  are  peculiar  to  two  dimen- 
sions. However,  the  Borel  sets  in  two  dimensions  can  be 
generated  from  the  diamond  sets,  and  possibly  a higher 


dimensional  development  could  proceed  from  a purely  set 
viewpoint  without  recourse  to  points.  Such  & development 
would  use  the  spectral  measure  E with  its  projections  E (A) 
and  not  even  try  to  bother  with  point  projections  such 
as  Ex.  We  feel  that  this  technique  would  end  up  produc- 
ing much  the  same  results  paralleling  classical  resolution 
space  as  in  the  two-dimensional  case. 

Although  most  of  the  results  mirrored  those  of 
classical  resolution  space,  one  new  type  of  operator  was 
obtained  which  appears  to  be  a generalization  of  the 
memoryless  operator.  This  is  the  spacelike  operator 
which  was  introduced  in  the  decomposition  theorem.  It 
is  conjectured  that  if  the  inverse  of  a spacelike  oper- 
ator exists,  then  it  will  also  be  spacelike. 

However,  other  than  this  one  interesting  conjecture 
the  extension  of  resolution  space  to  a relativistic  set- 
ting seems  to  offer  little  promise  of  new  insights  into 
operator  theoretic  systems  theory. 
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